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EXECUTIVE  SUMMARY 


This  report  describes  a  software  tool  for  two-way  analysis  of  variance  (ANOVA) 
with  repeated  measures  on  one  factor  for  use  in  a  multitude  of  problems,  including  the 
analysis  of  Affymetrix  GeneChip™  microarray  data.  The  proposed  software  has  been 
used  to  analyze  more  than  22,000  probe  sequences  in  less  than  1  minute. 

The  tool  is  entirely  written  in  Java  and,  as  a  result,  is  platform-independent.  The 
current  implementation  of  the  tool  only  allows  for  one-way  and  two-way  ANOVA. 
However,  the  internal  data  structure  is  designed  to  be  able  to  hold  data  for  multi-way 
ANOVA.  The  output  is  written  to  a  tab-delimited  text  file  and,  accordingly,  the  tool  can 
be  easily  extended  to  save  the  results  directly  into  a  relational  database. 
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INTRODUCTION 


The  development  of  increasingly  pot/verful  DNA  microarrays  for  the  study  of  gene 
expression  has  created  a  need  for  powerful  computational  tools  with  which  to  analyze 
the  data.  For  example,  a  late-generation  Affymetrix  oligonucleotide  array  (the  U1 33A) 
contains  probe  sets  for  22,283  sequences.  Each  of  these  is,  in  theory,  capable  of 
detecting  a  change  in  expression  for  the  sequence  it  detects  independently  of  all  other 
probe  sets  on  the  array.  This  means  that  any  statistical  test  designed  to  detect 
individual  sequence  changes  in  expression  will  need  to  be  repeated  22,283  times  over, 
and  the  output  will  need  to  be  formatted  in  a  manner  that  is  easily  imported  into 
software  packages  capable  of  handling  large  amounts  of  data,  such  as  relational 
databases. 

When  simple  experimental  designs  are  involved,  such  as  (for  example)  five  to  ten 
paired  sets  of  samples,  the  mathematical  analyses  are  fairly  straightforward  and  can  be 
performed  readily  with  a  commercially  available  spreadsheet,  such  as  Microsoft  Excel, 
which  can  perform  a  T-test  and  output  the  'esulting  P-value  into  the  same  row  as  the 
sequence  under  from  which  it  was  generated.  For  more  complex  research  designs, 
however,  such  as  comparisons  between  three  or  more  means  or  two-factorial  designs 
(e.g.,  multiple  experimental  conditions  tested  at  varying  time  points),  the  use  of 
spreadsheets  becomes  increasingly  impractical  due  both  to  computational  constraints 
(memory,  processing  power)  and  the  difficulty  of  managing  the  output.  For  example,  in 
Microsoft  Excel,  a  one-way  analysis  of  variance  (ANOVA)  generates  a  table  of  output, 
which  greatly  complicates  the  process  of  maintaining  a  one-to-one  relationship  between 
the  input  data  and  the  output  of  the  calculation. 

A  need  exists  for  computational  tools  that  can  rapidly  perform  ANOVA  on  DNA 
microarray  data.  ANOVA  is  an  excellent  method  for  analyzing  microarray  data  when 
multiple  groups  or  factors  are  involved  (such  as  treatment  vs.  time),  as  it  permits 
statistical  evaluation  of  changes  in  gene  expression  and  variability  in  the  array  data 
(3:4).  ANOVA  has  been  productively  used  to  calculate  the  contribution  and  significance 
of  different  treatments  to  changes  in  gene  expression  in  published  literature  (e.g.,  see 
reference  (1)), 

There  are  several  software  tools  that  can  perform  ANOVA,  including  commercial 
programs,  such  as  SAS  (6)  and  SPSS  (8),  as  well  as  free  tools,  such  as  the  R 
programming  language  (1 1 ).  However,  the  vast  amount  of  data  generated  by  even  a 
comparatively  small  microarray  experimeni  can  unmask  software  limitations  that  are  not 
always  apparent  when  the  number  of  variacles  under  consideration  is  small.  For 
example,  SPSS  generates  large  amounts  of  output  other  than  P  values  for  each 
ANOVA  calculation:  this  feature  is  very  useful  when  the  number  of  independent 
calculations  is  small,  but  cumbersome  to  deal  with  when  attempting  to  perform  the 
same  computation  on  hundreds  to  thousands  of  separate  variables.  Furthermore,  the 
Windows  Graphical  User  Interface  versions  of  SPSS  available  when  this  project  was 
started  were  not  capable  of  performing  thousands  of  two-way  ANOVA  calculations  at 
one  sitting  on  the  hardware  platforms  available  at  the  time.  Conversely,  console  / 
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character-based  statistical  programming  languages,  such  as  R  and  SAS,  are  highly 
capable  of  performing  the  required  analyses  but  have  their  own  unique  disadvantages. 
For  example,  the  cost  of  a  license  to  run  a  general-purpose  program,  such  as  SAS, 
makes  it  an  unattractive  option  for  a  highly  specialized  requirement.  The  R 
programming  language  is  available  as  freeware,  but  also  suffers  the  performance 
limitations  inherent  to  a  general-purpose,  interpreted  language,  which  greatly  limits  the 
speed  at  which  repetitive  calculations  can  be  performed.  SPSS  also  has  a  built-in 
scripting  language,  but  unlike  R,  which  is  available  as  freeware,  a  script  written  in  SPSS 
is  of  value  only  to  users  who  purchase  a  license  to  run  SPSS. 

We  therefore  decided  to  develop  and  implement  a  software  tool  that  allows  the 
user  to  perform  one-way  and  two-way  ANOVA  with  or  without  repeated  measures  for 
any  number  of  instances.  We  chose  to  write  the  program  in  Java  to  maximize  the 
number  of  platforms  on  which  it  can  execute.  The  tool  generates  results  as  tab- 
delimited  files  that  can  be  readily  imported  into  relational  databases. 

This  software  tool  will  be  maintained  by  the  Bioinformatics  Cell  of  the  U.S.  Army 
Medical  Research  and  Materiel  Command,  Fort  Detrick,  Maryland. 

SAMPLE  DATA  SET 

Although  the  software  allows  the  user  to  carry  out  two-way  repeated  measures 
ANOVA  for  any  data  set,  it  was  written  with  a  particular  experiment  in  mind.  The 
experimental  design  has  three  conditions  (heat,  cold,  control)  in  which  gene  expression 
signals  are  measured  at  time  points  T  =  0  and  T  =  3.  The  experiment  is  repeated  six 
times,  and  the  data  from  each  trial  are  numbered  from  1  to  6.  There  are  a  total  of 
22,283  gene  sequences  in  each  data  set,  whose  signals  are  measured  in  each 
experiment.  Each  gene  sequence  is  labeled  in  the  database  with  a  unique  identifier 
(ID).  An  abstract  representation  of  this  data  set  for  a  single  gene  is  shown  in  Table  1 
(next  page). 
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T able  1 :  Abstract  representation  of  a  sample  data  for  a  single  gene,  where  Yij^  is  the 
measured  expression  of  a  gene  in  sample  /  for  condition  j  at  time  k.  Each  gene  is 
measured  in  6  samples  for  3  conditions  at  2  time  points. 


Condition 

Sample 

Time 

T=0 

T=3 

Heat 

1 

Yu, 

Y11? 

2 

y2ii 

Ypi? 

3 

Y3U 

Y312 

4 

Y411 

Y412 

5 

Y511 

Ym? 

6 

Yen 

Y612 

Cold 

1 

Y121 

Yi?? 

2 

y2?i 

Y22? 

3 

y32i 

Y.122 

4 

Ya?, 

Y4?? 

5 

y52i 

Y522 

6 

y62i 

Vfi.s? 

Control 

1 

Yi31 

Yi32 

2 

Y231 

Y?32 

3 

Y331 

Y332 

4 

y43i 

Y43? 

5 

Y531 

Y53P 

6 

Y63, 

Yes? 
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MATERIALS  AND  METHODS 

TWO-WAY  ANOVA  WITH  REPEATED  MEASURES  ON  ONE  FACTOR 

The  goal  of  our  data  analysis  is  to  perform  two-way,  repeated  measures  ANOVA 
for  each  gene  sequence.  In  our  sample  dataset,  the  between-subjects  factor  is 
condition,  and  the  within-subjects  repeated  measures  factor  is  time.  For  each  gene 
sequence  in  the  dataset,  the  following  statistics  are  calculated: 

■  P-value  for  the  main  effect  {condition). 

■  P-value  for  the  within-subjects  factor  (time). 

•  P-value  for  the  interaction  {time  *  condition). 

ONE-WAY  ANOVA 

In  addition  to  the  two-way  repeated  measures  ANOVA,  an  additional  secondary 
analysis  is  also  performed.  Because  one  of  the  factors  is  repeated,  it  is  possible  to 
derive  a  slope  of  the  factor.  Each  gene  is  measured  at  two  different  time  points  in  three 
different  conditions,  which  allows  us  to  calculate  a  slope  for  each  gene  in  each 
condition. 

■  Slope{heat)  =  [{signal(heat)  at  T=3)  -  (signal(heat)  at  T=0)]  /  [3  -  0] 

■  Slope(cold)  =  [(signal(cold)  at  T=3)-  (signal(cold)  at  T=0)]  /  [3  -  0] 

■  Slope(control)  =  [(signal(control)  at  T=3)  -  (signal(control)  at  T=0)]  /  [3  -  0] 

In  other  words,  for  each  gene  in  each  of  the  six  experiments,  three  derived 
variables  are  computable:  Slope(heat),  Slope(cold),  Slope(control).  Using  these  slope 
values  we  can  perform  a  one-way  ANOVA  with  condition  (heat,  cold,  control)  as  the 
between-groups  factor  for  each  gene  in  the  dataset.  We  compute: 

■  P-value  for  the  main  effect  {condition). 

■  Post-hoc  Tukey's  test  for  the  comparisons: 

■  Heat  vs.  Control 

•  Heat  vs.  Cold 

■  Cold  vs.  Control 

The  following  sections  describe  in  detail  the  one-way  and  the  two-way,  repeated 
measures  ANOVA  models  used  for  computing  these  results. 

TWO-WAY  ANOVA  MODEL  WITH  REPEATED  MEASURES  ON  ONE  FACTOR 

Suppose  two  factors,  A  and  B,  are  involved  in  a  study  and  repeated  measures 
are  taken  on  B,  where  A  has  a  levels,  B  has  b  levels,  and  n  subjects  are  available  for 
each  level  of  A.  Let  Yijk  denote  the  observation  for  the  i-th  subject  when  A  is  at  level  j 
and  B  at  level  k.  Let  a,  denote  the  main  effect  of  A  at  level  j,  Pi<  the  main  effect  of  B  at 
level  k,  and  (a^)jkthe  interaction  effect,  when  A  is  at  level  j  and  B  at  level  k,  where,  j=1, 
2,  ...,  a;  and,  k=1,  2,  ....  b.  In  addition,  let  p,(j^  denote  the  subject’s  main  effect. 
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Assuming  that  no  interactions  exist  between  treatments  and  subjects,  we  have  the 
following  model  (5): 


=  A...  +  Ao)  +£ijk  (1) 

where,  i)  [j...  is  a  constant:  ii)  piQ)  is  iid  N(0,Op^)  iii)  O/,  ^k,  and  (ap)jk  are  constants,  such 
that  I  Oj  =  d,lBk  =  0,  Ij(a(S)jk  =  0  for  all  k,  and  Tk(ctP)jk  =  0  for  all  j:  iv)  the  error  terms  % 
are  iid  N(0,  cr);  and  v)  p,(j;  and  e^k  are  independent. 


In  our  example  dataset,  for  each  gene,  we  need  to  perform  the  two-way  repeated 
measures  ANOVA,  where  the  between-subjects  factor  is  the  condition  (heat,  cold  or 
control)  and  the  within-subjects  repeated  measures  factor  is  time. 

The  effect  of  each  factor  and  their  interactions  are  calculated  as  (5): 


MSA= - ^ - — 

(2) 

MSB= - ^ - 

b-\ 

(3) 

(4) 

- 

a{n-\) 

(5) 

msb.s(a)  =  '  - 

(6) 

Iizn* 

Y  -  '  '  ' 

^  -  nab 

(7) 

^  “iicT  =  independently  and  identically  distributed 
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(8) 


II7, 


tjk 


y  =  '  * 

ZIK 


ijk 


Y..r 


na 


7. 


•-'*  n 


(9) 


(10) 


(11) 


where  /=  7,  2,  ...,  n;  j  =1,2,  ...,  a;  and  k  =  1,2,  6;  n  is  the  number  of  times  the 

experiment  is  performed,  a  is  the  number  of  different  values  the  factor  for  the  main 
effect  can  have,  and  b  is  the  number  of  different  values  for  the  repeated  factor.  In  the 
case  of  our  example  data  set,  n  =  6,  a  =  3,  b  =  2,  and  Yijk  is  the  expression  value  of  the 
gene  in  i-th  sample  in  j-th  condition  at  k-th  time  point. 


The  p-value  for  the  main  effect  (i.e.,  the  factor  condition  in  our  example)  is 
calculated  as  the  probability  of  F  >F*,  where  F  is  a  random  variable  with  numerator 
degrees  of  freedom  equal  to  a  -  f,  denominator  degrees  of  freedom  equal  to  a(n-  1), 
and 


F* 


MSA 

mss{a)  ■ 


(12) 


The  p-value  for  the  within-subjects  factor  (i.e.,  the  factor  time  in  our  example)  is 
calculated  as  the  probability  of  F  >F*,  where  F  is  a  random  variable  with  numerator 
degrees  of  freedom  equal  to  b-  1,  denominator  degrees  of  freedom  equal  to  a(n  -  1)(b 
-  1),  and 


MSB 

MSB.S{a)  ■ 


(13) 


The  p-value  for  the  interaction  is  calculated  as  the  probability  of  F  >F*,  where  F 
is  a  random  variable  with  numerator  degrees  of  freedom  equal  to  (a  -  1)(b  -  1), 
denominator  degrees  of  freedom  equal  to  a(n  -  1)(b  - 1),  and 
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MSAB 


(14) 


MSB.S{a)  ■ 

ONE-WAY  ANOVA  MODEL 


Let  us  suppose  there  is  only  one  fac:or  involved  in  a  study,  where  the  factor  has  r 
levels,  and  n,  denotes  the  number  of  cases  for  the  i-th  factor  level.  The  single-factor 
ANOVA  model  can  be  stated  as  follows  (5): 

Yy=u^+e,^  (15) 

where  (i)  Y/j  is  the  observation  in  the  j-th  trial  for  the  i-th  factor  level  (for  i=  1,  2,  r  and 
j  =1,2,...,  nA;  (ii)  represents  the  mean  response  for  the  i-th  factor  level;  and  (iii)  £,y 
are  iid  N(0,  cr). 

As  discussed  before,  in  the  example  dataset,  for  each  gene  in  each  condition,  a 
signal  slopes  can  be  computed  and  used  for  the  one-way  ANOVA  computation  with  the 
condition  as  the  between-groups  factor.  The  signal  slope  for  each  gene  in  a  given 
condition  and  a  given  sample  can  be  calculated  as: 


T  -T 
3-0 


(16) 


where  the  signal  in  condition  i  in  sample  /at  time  =  3,  and  is  the  signal  in 
condition  /  in  sample  /  at  time  =  0. 

After  calculating  the  slopes,  the  data  for  each  gene  in  Table  1  can  be 
represented  as  shown  in  Table  2. 


Table  2:  Representation  of  the  data  for  a  gene  after  calculating  the  signal  slopes  for  the 
given  condition  and  the  given  sample,  where  Yij  is  the  slope  for  condition  /  and  sample  /. 
For  the  example  data  set,  we  can  calculate  six  slopes  for  each  of  the  three  conditions. 


Condition 

Slopes 

Heat  (ui) 

Yi2 

Yi3 

Yu 

^75 

!  Yre 

Cold  (U2) 

y2i  : 

Y,7 

'^23 

V24 

V25 

Y76 

Control  (u:^) 

Y31 

Y32 

Y33 

Y34 

^35 

Yse 

Now  we  can  calculate  the  p-values  for  the  main  effect  (i.e.,  the  condition  in  our 
example)  using  Equations  17  to  22  (5),  defined  as  follows; 


MSTR  =  ^  ^  ’ 

r-1 


(17) 
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MSE  = 


(18) 


TL^.rY)' 

rij-r 


1 7. 


7,-, 


7 


y=i 


n. 


ZZ7, 


«r  =Z"' 


(19) 

(20) 
(21) 


where  /  =  1 ,  2 . r;  y  =  1 ,  2 . nf,r  Is  the  number  of  different  values  the  factor  for  the 

main  effect  can  have,  and  n,  is  the  maximum  number  of  slopes  for  each  different  value 
of  the  factor  for  the  main  effect.  For  our  example  data  set,  r  =  3  and  ni  =  n2  =  03  =  6. 
Equation  19  calculates  the  average  slope  for  each  gene  in  each  condition,  and  Equation 
20  calculates  the  overall  average  slope  for  each  gene. 

Using  the  F-distribution,  the  p-value  for  the  main  effect  (i.e.,  the  factor  condition 
in  our  example)  is  calculated  as  the  probability  of  F  >F*,  where  F  is  a  random  variable 
with  numerator  degrees  of  freedom  equal  to  (r-  1),  denominator  degrees  of  freedom 
equal  to  nj-r,  and 


_  MSTR 
~  MSE  ■ 


(22) 


TUKEY’S  HSD  COMPARISONS 

The  second  goal  of  the  one-way  ANOVA  is  to  do  Post-hoc  Tukey's  HSD 
(Honestly  Significant  Difference)  comparisons  for  all  pairs.  However,  instead  of 
calculating  the  p-value  for  each  pair,  the  program  provides  confidence  limits  with  family 
confidence  coefficient  of  at  least  1  -  a.  The  confidence  limits  provide  estimation  of  the 
range  over  which  the  difference  of  the  mean  value  for  the  pair  varies  with  the  family 
confidence  coefficient  of  at  least  1  -  a.  If  the  confidence  limits  contain  zero,  the  means 
of  the  two  conditions  are  the  same  at  some  point.  However,  if  the  limits  does  not  include 
zero,  the  means  are  never  the  same  with  the  p-value  of  a.  For  our  example  data, 
consider  that  for  a  given  gene,  the  mean  signal  value  in  condition  heat  is  p?  and  in 
condition  cold  is  \J2.  Assume  that  the  confidence  limits  (calculated  as  /l/j  -  /J2)  are  [-100, 
100]  and  a  =  0.05.  Since,  the  limits  include  zero,  we  can  say  that  the  means  ^J1  and  ^J2 
are  not  different  with  95%  confidence.  However,  if  the  confidence  limits  were  either  [- 
100,  -50]  or  [100, 1000],  we  can  say  that  the  means  jJi  and  iJ2  are  different  with  95% 
confidence  (p-value  of  0.05). 
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The  confidence  limits  for  ali  pair-wise  comparisons  D  =  jU?  -  /1/2  with  family 
confidence  coefficient  of  at  least  1  -  a  are  d+t-sid},  where 


D=V  - 


YrY> 


5{Z)}  = 


T  =  -j=q{\-a-,r,nj-r) 

■Jl 


(23) 

(24) 

(25) 


In  the  case  of  our  example  dataset,  r=  3,  nr  =  18,  and  a  of  0.01 , 0.05,  and  0.1, 
the  function  q(1  -  a;  r,  nj-r)  in  Equation  25  evaluates  to  4.83,  3.67,  and  3.14, 
respectively. 

PROGRAM  DESCRIPTION 


The  program  is  written  in  Java.  As  such,  it  is  platform  independent  and  runs  on 
Windows,  Linux,  Mac  OS,  and  Unix.  To  use  the  tool,  the  following  files  are  required: 

The  ANOVA  Java  program  file  itself,  compiled  as  an  executable  Java  program;  a 
configuration  file  (described  below,  with  an  example  shown  in  Appendix  A)  that  supplies 
the  program  with  the  basic  information  it  needs  to  perform  the  computation  and  allows 
the  user  to  specify  various  parameters  and  the  fonnat  of  the  input  file;  the  data  file  to  be 
analyzed;  and  a  mathematical  library  publicly  available  by  the  Jakarta  project  (10), 
which  is  used  to  calculate  the  F  distribution. 


For  users  experienced  in  compiling  Java  programs,  the  source  code  is  listed  in 
Appendix  B  to  this  report.  Executable  code  with  sample  input  and  configuration  files 
can  also  be  obtained  directly  from  the  authors  (E-mail;  iaques.reifman@us.armv.mil. 
larrv  sonna@vahoo.com  or  purvesh@cs.wavne.edu  ).  A  minimally-annotated 
configuration  file,  test  data,  and  output  files  that  can  be  used  to  verify  program 
execution  are  included  in  Appendix  C. 

How  to  Use  the  Tool 


In  order  to  use  the  tool,  Java  must  be  installed  on  the  user’s  computer.  Java  can 
be  downloaded  from  Sun  Microsystems  (9).  The  following  command  can  be  used  at  the 
command-line  to  run  the  tool: 

java  anova.client.AnovaAnalysisI  <configuration  file  name>  -  d  <#  ofdigits> 

The  switch  “-d”  is  optional  and  is  followed  by  the  number  of  digits  that  the  user 
wishes  the  program  to  compute  for  each  P  value.  The  default  is  8-digit  precision. 

Note  that  the  Jakarta  commons-math  library  should  be  in  the  “classpath”  (the 
directory)  so  that  Java  can  find  it  during  execution.  The  classpath  can  be  specified  at 
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execution  time  by  “-cp”  switch.  For  instance,  if  commons-math  library  and  the  ANOVA 
code  is  in  the  same  directory,  the  following  command  can  be  used: 

java  -cp  “Jcommons-math-l .  1  Jar;./" anoava.client.AnovaAnalysisI 
<configuration  file  name>  (for  Windows) 


or 

java  -cp  “./commons-math-l .1  .jar:./’  anova.client.AnovaAnalysisI 

<configuration_fi!e__name>  (for  Unix/Linux) 

To  start  a  command-line  program  on  a  Microsoft  windows  platfomn,  click  on 
“Start,”  and  then  “Run.”  Type  “cmd”  in  the  dialog  box  that  pops  up  and  press  the  “OK” 
button.  A  DOS  window  will  open;  the  user  will  then  need  to  change  the  directory  (using 
the  “cd”  command)  to  the  directory  where  the  ANOVA  tool  is  located.  Once  in  that 
directory,  the  commands  listed  above  can  be  entered  to  start  the  program.  For  the 
convenience  of  inexperienced  Windows  users,  a  batch  program  that  executes  the 
ANOVA  tool  is  included  in  the  electronic  distribution  (available  from  the  authors). 

On  a  Unix  or  Linux  workstation,  right-clicking  the  mouse  and  selecting  “terminal” 
or  “console”  from  the  pop-up  menu  usually  opens  a  command-line  terminal. 

The  configuration  file  should  be  prepared  as  described  in  the  following  section. 
The  results  are  saved  in  files  as  specified  in  the  configuration  file.  The  formats  of  the 
input  file  and  the  output  file  are  described  below. 

Configuration  File 

Use  of  the  ANOVA  Java  tool  requires  that  the  user  supply  the  program  with  a 
configuration  file.  This  file  gives  the  program  the  basic  information  it  needs  to  properly 
perform  the  desired  computation.  A  sample  file  is  included  in  Appendix  A,  and  the 
required  contents  of  the  file  are  listed  in  Tables  3  and  4  (below).  The  ANOVA  Java  tool 
requires  that  the  configuration  file  be  a  simple  ASCII  file;  thus,  if  the  configuration  file  is 
created  using  Microsoft  Word  or  Notepad,  it  must  be  saved  as  a  plain  “text”  file. 

The  format  of  the  configuration  file  is  similar  to  that  of  windows  “.ini”  files.  The 
configuration  file  is  divided  into  various  sections,  called  “paragraphs.”  Each  paragraph 
starts  with  the  name  of  the  paragraph  enclosed  between  brackets,  “[“  and  “]”.  This 
approach  allows  for  easy  addition  of  the  new  paragraphs  in  the  future  while  keeping  the 
new  configuration  file  compatible  with  its  previous  versions. 

Any  line  starting  with  ”#”  in  the  configuration  file  is  treated  as  a  comment  line,  i.e., 
does  not  contain  instructions  that  the  program  will  use. 

The  configuration  parameters  are  divided  into  two  parts:  (1)  global  configuration 
parameters  and  (2)  factor  configuration  parameters.  The  global  configuration 
parameters  are  specified  in  the  “config”  paragraph,  and  each  ANOVA  factor  is 
described  in  a  “Factor/V”  paragraph,  where  N  is  an  integer.  Note  that  a  “Factor” 
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paragraph  refers  to  each  index  used  to  describe  a  data  item,  not  the  factors  (in  a 
statistical  sense)  used  in  the  analysis.  For  instance,  in  the  case  of  the  sample  data  set 
in  Table  1,  each  data  Item  is  described  by  three  indices:  (1)  a  sample  number  (1-to-6), 
(2)  a  condition  (heat,  cold,  or  control),  and  (3)  a  time  point  (t=0  or  t=3).  Of  these  three 
indices,  only  condition  and  time  indices  refer  to  the  two  ANOVA  factors  used  in  the 
analysis.  However,  in  order  for  the  program  to  know  which  data  to  pair  up  in  time,  the 
sample  number  factor  is  required.  The  integer  number  used  to  describe  each  “Factor” 
paragraph  must  start  at  1  and  should  be  incremented  by  1  for  each  subsequent  index. 
For  instance,  in  the  case  of  the  sample  data  set,  there  are  three  indices  used  to 
describe  a  data  item.  Therefore,  the  configuration  file  must  contain  three  “Factor” 
paragraphs,  i.e.,  Factorl,  Factore2,  and  Factors. 

Each  paragraph  (either  “config”  or  “FactorN”)  assigns  value  to  one  or  more 
parameters.  The  syntax  to  assign  a  value  to  a  parameter  is  “name  =  value."  The 
parameters  for  a  paragraph  can  be  specified  in  any  order.  The  paragraph  name  and  the 
parameter  name  are  case-insensitive. 

The  “config”  paragraph  contains  the  oarameters  described  in  Table  3,  and  the 
“Factor”  paragraph  parameters  are  described  in  Table  4.  Appendix  A  provides  a  sample 
configuration  file. 


Table  3.  Description  of  the  parameters  for  the  “config”  paragraph. 


Parameter  Name 

Description 

NumOfFactors 

The  total  number  of  indices  used  to 
describe  each  data  item.  It  must  be  an 
integer  number.  For  instance,  in  the  case 
of  our  sample  data  set.  NumOfFactors  =  3. 

DataFile 

The  location  of  the  input  file  containing  the 
data.  The  location  can  be  either  an 
absolute  path  or  a  relative  path. 

TwoWayOutputFile 

File  name  in  which  the  results  of  the  two- 
way  ANOVA  should  be  stored. 

HeadersForOneWayAnova 

Factor  name  that  should  be  used  for 
column  headers  in  the  file  specified  by 
OneWayOutputFile  (see  below). 

OneWayOutputFile 

File  name  in  which  the  results  of  the  one¬ 
way  ANOVA  should  be  stored. 
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Table  4:  Description  of  parameters  for  the  “factor”  paragraph. 


Parameter  Name 

Description 

Name 

Name  of  the  index  described  by  the 
paragraph.  It  can  be  any  free-form  string. 

Column 

The  column  number  in  the  input  file  in 
which  the  value  for  the  index  can  be  found. 
The  first  column  in  the  input  file  has  index 

1 .  The  input  file  is  the  file  specified  by 
DataFiie  parameter  in  the  "config" 
paragraph  (See  Table  3). 

MaxNumOfIndex 

The  maximum  value  possible  for  the  index. 

Enumerate 

Optional.  If  an  index  has  categorical 
values,  each  category  can  be  assigned  an 
integer,  starting  at  1  and  incremented  by  1 . 
When  enumerating  an  index,  the  category 
name  and  its  index  must  be  separated  by 
a  tab  character.  See  Appendix  A  for  an 
example. 
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RESULTS 


INPUT  FORMAT 

The  input  file  is  specified  via  the  DataFile  parameter  of  the  “config”  paragraph 
(see  Table  3).  The  data  for  ANOVA  is  always  a  matrix,  where  each  element  of  the 
matrix  is  a  number.  Therefore,  the  format  of  the  input  file  for  any  ANOVA  software 
requires  that  the  input  file  describes  one  matrix  element  per  line.  For  an  example  of  the 
input  file  format,  see  Appendix  A. 

OUTPUT  FORMAT 

The  ANOVA  tool  carries  out  one-way  and  two-way  ANOVA  computations  on  the 
data  in  the  input  file.  The  names  of  the  output  files  for  both  of  the  analyses  are 
specified  by  OneWayOutputFile  and  TwoWayOutputFile  parameters  of  the  “config” 
paragraph  (see  Table  3).  The  formats  of  the  one-way  ANOVA  output  file  and  the  two- 
way  ANOVA  output  file  are  described  in  Appendix  A. 

DESCRIPTION  OF  APPENDIX  B 

Appendix  B  provides  a  listing  of  the  Java  program  along  with  its  related  routines, 
a  list  of  the  configuration  file  for  the  sample  data  set  discussed  in  the  report,  and  an 
input  file  for  20  genes  with  its  associated  one-way  and  two-way  ANOVA  output  files. 
Appendix  C  contains  a  test  dataset  and  the  output  resulting  from  application  of  the 
ANOVA  tool  to  this  dataset. 


COMPARISON  OF  OUTPUT  GENERATED  BY  THE  JAVA  PROGRAM  TO  THE 
OUTPUT  GENERATED  BY  SPSS  10.0.5 

To  validate  the  Java  software  routine  described  in  this  report,  we  generated  a 
test  file  of  90  sequences  from  our  sample  dataset  and  computed  the  two-way,  repeated- 
measures  ANOVA  P  values  for  condition,  time,  and  the  interactive  factor  condition*time 
using  SPSS  1 0.0.5.  The  output  was  rounded  to  1 5  significant  digits  (1 6  digits  total, 
including  the  leading  zero).  The  sample  file  used  to  validate  the  output  against  SPSS 
contained  different  data  from  the  file  used  to  code  and  debug  the  software  tool 
(Appendix  B).  To  insure  that  the  computational  algorithm  would  hold  true  not  only  for 
mid-range  P  values  but  also  for  extremes,  we  included  genes  in  this  validation  test  file 
that  had  P  values,  as  measured  by  the  Java  program,  both  very  close  to  zero,  as  well 
as  close  to  1 .  To  verify  that  results  would  be  comparable  across  hardware  platforms, 
the  computation  with  the  Java  program  was  performed  on  a  different  computer  than  the 
one  on  which  the  SPSS  computation  was  performed. 

The  results  are  shown  in  Table  5.  The  differences  between  the  two  computations 
were  negligibly  small  and  within  the  margin  of  error  expected  from  floating-point 
computation  error  (2).  The  comparison  showed  that  the  Java  program  output  is 
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essentially  the  same  as  that  obtained  with  one  of  the  major  commercial  software 
statistical  packages. 


Table  5.  Comparison  of  results  generated  by  the  JAVA  ANOVA  program  to  the  results 
generated  by  SPSS  10.0.5  No  more  than  four  significant  digits  are  reported  in  this 

tabie. 

Table  5A.  Raw  data 


JAVA  PROGRAM 

RVALUES 

SPSS  10.0.5 

RVALUES 

DIFFERENCE 

Probe  ID 

COND 

TIME 

TIME 

*COND 

I  COND 

TIME 

TIME 

*COND 

,  COND 

TIME 

TIME  ‘  COND 

200878  at 

0.0362 

0.0859 

0.6661 

0,0362 

0.0859 

0.6661 

4.305E-11 

-1.530E-11  -6.542E-12 

201060  X  at 

0.2106 

0.3122 

0.0329 

0.2106 

0.3122 

0.0329 

-6.273E-12 

-1.475E-11  6.751E-11 

201118  at 

0.3434 

0.0619 

0.1243 

0.3434 

0.0619 

0.1243 

5.129E-12 

-1.426E-11  -3.784E-11 

201162  at 

0.6202 

0.5742 

0.0001 

0.6202 

0.5742 

0.0001 

-2.817E-11 

2.968E-12  1.683E-12 

201163  s  at 

0.4764 

0.4796 

0.6164 

0.4764 

0.4796 

0.6164 

-9.416E-11 

9.589E-11  -3.160E-11 

201206  s  at 

0.7017 

0.6065 

0.6438 

0.7017 

0.6065 

0.6438 

-1.904E-12 

8.835E-13  -1.352E-11 

201858  s  at 

0.3759 

0.3434 

0.4293 

0.3759 

0.3434 

0.4293 

2.426E-11 

-4.325E-12  4.860E-12 

201909  at 

0.2241 

0.7017 

0.8814 

0.2241 

0.7017 

0.8814 

2.836E-10 

1.355E-12  3.899E-12 

202018  s  at 

0.4326 

0.6594 

0.9795 

0.4326 

0.6594 

0.9795 

4.400E-12 

7.755E-12  -1.266E-14 

202611  s  at 

0.1436 

0.7317 

0.0020 

0.1436 

0.7317 

0.0020 

-1.394E-11 

2.370E-13  3.665E-12 

203153  at 

0.5695 

0.0174 

0.2065 

0.5695 

0.0174 

0.2065 

-6.065E-12 

1.328E-11  -7.408E-12 

203290  at 

0.6791 

0.1096 

0.5627 

0.6791 

0.1096 

0.5627 

-4.216E-12 

2.089E-10  -7.493E-12 

203828  S  at 

0.1447 

0.2542 

0.0073 

0.1447 

0.2542 

0.0073 

-1.320E-11 

-1.450E-10  7.303E-11 

203979  at 

0.0826 

0.0462 

0.1427 

0.0826 

0.0462 

0.1427 

3.028E-12 

-8.394E-11  -1.456E-11 

204174  at 

0.6158 

0.9531 

0.3105 

0.6158 

0.9531 

0.3105 

-3.216E-11 

2.551  E-1 3  1.500E-11 

204351  at 

0.6527 

0.2599 

0.9837 

0.6527 

0.2599 

0.9837 

-1.016E-11 

-1.156E-10  -2.442E-15 

204439  at 

0.6672 

0.1829 

0.4941 

0.6672 

0.1829 

0.4941 

-6.307E-12 

4.619E-11  -5.690E-11 

204475  at 

0.9726 

0.5315 

0.3444  1 

0.9726 

0.5315 

0.3444 

-8.116E-14 

1.542E-11  4.962E-12 

204669  s  at 

0.9832 

0.7536 

0.8736 

0.9832 

0.7536 

0.8736 

-3.886E-15 

9.881E-14  6.278E-12 

204750  s  at 

0.0173 

0.7846 

0.1477 

0.0173 

0.7846 

0.1477 

2.240E-12 

-7.905E-12  -1.135E-11 

204939  s  at 

0.2446 

0.1787 

0.5278 

0.2446 

0.1787 

0.5278 

1.382E-10 

5.678E-11  -2.143E-11 

204970  S  at 

0.5543 

0.0079 

0.6932 

0.5543 

0.0079 

0.6932 

-9.683E-12 

3.952E-12  -2.577E-12 

205000  at 

0.1022 

0.1543 

0.0019 

0.1022 

0.1543 

0.0019 

-1.331E-10 

1.442E-11  4.391E-12 

205033  S  at 

0.7948 

0.2621 

0.1250 

0.7948 

0.2621 

0.1250 

6.251E-12 

-1.061E-10  -3.653E-11 

205118  at 

0.9538 

0.7522 

0.7113 

0.9538 

0.7522 

0.7113 

-1.966E-12 

8.749E-14  -1.334E-12 

205254  X  at 

0.2553 

0.5994 

0.5408 

0.2553 

0.5994 

0.5408 

9.550E-11 

1.118E-12  -1.456E-11 

205403  at 

0.3065 

0.0158 

0.3602 

0.3065 

0.0158 

0.3602 

1.713E-11 

2.152E-11  3.885E-11 

205557  at 

0.4041 

0.4963 

0.0513 

0.4041 

0.4963 

0.0513 

1.041  E-11 

5.390E-11  7.260E-12 

205576  at 

0.3173 

0.0357 

0.0098 

0.3173 

0.0357 

0.0098 

1.201E-11 

2.171E-12  2.371E-11 

206385  s  at 

0.1824 

0.2003 

0.6906 

0.1824 

0.2003 

0.6906 

-1.990E-11 

2.018E-11  -2.827E-12 

206522  at 

0.5928 

0.0976 

0.4570 

0.5928 

0.0976 

0.4570 

-2.908E-12 

-6.144E-12  -1.619E-10 
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206582 

s  at 

0.7456 

0.3568 

0.0095 

0.7456 

0.3568 

0.0095 

4.432E-11  -2.570E-12  2,668E-11 

206676 

at 

0.2883 

0.7038 

0.8468 

0.2883 

0.7038 

0.8468 

3.130E-11  1.128E-12  4.777E-13 

206697 

s  at| 

0.0257 

0.3305 

0.5646 

0.0257 

0.3305 

0.5646 

1.824E-12  -7.179E-12  -7.059E-12 

206700 

s  at 

0.0000 

0.0311 

0.0225 

0.0000 

0.0311 

0.0225 

-8.912E-12  4.630E-12  3.372E-12 

206785 

s  at 

0.9406 

0.0336 

0.0010 

0.9406 

0.0336 

0.0010 

2,520E-14  3.014E-12  1.086E-11 

206871 

at 

0.3615 

0.0539 

0.3966 

0.3615 

0.0539 

0.3966 

3.737E-11  -3.370E-11  1.302E-11 

206877 

at 

0.0494 

0.2246 

0.0004 

0.0494 

0.2246 

0.0004 

8.887E-12  6.588E-12  6.510E-12 

207177 

at 

0.6558 

0.3032 

0.6550 

0.6558 

0.3032 

0.6550 

-9,170E-12  -2,096E-11  -9.434E-12 

207269 

at 

0.7231 

0.8016 

0.7690 

0.7231 

0.8016 

0,7690 

-8.561E-13  -3.171E-12  1.824E-11 

207384 

at 

0.6949 

0.1093 

0.6136 

0.6949 

0.1093 

0.6136 

-2.429E-12  2.130E-10  -3.435E-11 

207521 

s  at 

0.3442 

0.7907 

0.7873 

0.3442 

0.7907 

0,7873 

4.994E-12  -5.646E-12  8.638E-12 

207802 

at 

0.3932 

0.6210 

0.4274 

0.3932 

0.6210 

0.4274 

1.444E-11  4.499E-13  5.144E-12 

207824 

s  at 

0.6045 

0.2386 

0.9744 

0.6045 

0.2386 

0.9744 

-1.985E-12  -2.721E-10  -5.351E-14 

207892 

at 

0.3600 

0.9718 

0.5086 

0.3600 

0.9718 

0.5086 

3.917E-11  -3,553E-14  -3.753E-11 

207935 

s  at 

0.1461 

0.4933 

0.6721 

0.1461 

0.4933 

0.6721 

-1.226E-11  5.980E-11  -5.356E-12 

208450 

at 

0.5385 

0.1307 

0.0387 

0.5385 

0.1307 

0.0387 

-1.556E-11  5.641  E-11  3.11  IE-11 

208470 

s  at 

0.0064 

0.2684 

0.5642 

0.0064 

0.2684 

0.5642 

3.905E-12  -8.262E-11  -7.150E-12 

208602 

X  at 

0.4534 

0.2518 

0,6643 

0.4534 

0.2518 

0.6643 

-1.788E-10  -1.598E-10  -6.938E-12 

208791 

at 

0.4954 

0.0331 

0,5268 

0.4954 

0.0331 

0.5268  ; 

-5.485E-11  3.291E-12  -2.207E-11 

209156 

s  at 

0.0112 

0.0123 

0,0528 

0.0112 

0.0123 

0.0528 

1.410E-11  6.701E-11  6.228E-12 

209395 

at 

0.0415 

0,7655 

0.9683 

0.0415 

0.7655 

0.9683 

2.193E-11  3,941E-14  -1.980E-13 

209602 

s  at 

0.9892 

0.2026 

0.4324 

0.9892 

0.2026 

0.4324 

O.OOOE+00  1.813E-11  4.417E-12 

209840 

s  at 

0.1848 

0.1449 

0,1737 

0.1848 

0.M49 

0.1737 

-1.800E-11  2.461  E-11  -2.883E-11 

210116 

at 

0.1871 

0.1538 

0.4936 

0.1871 

0.'538 

0.4936 

-1.634E-11  1,482E-11  -5.774E-11 

210321 

at 

0.5831 

0,0343 

0.7447 

0.5831 

0.0343 

0.7447 

-3.957E-12  2,695E-12  4.569E-11 

210362  X  at 

0.5903 

0.3790 

0.3951 

0.5903 

0.3790 

0.3951 

-3.152E-12  -a.IME-ll  1.365E-11 

210690 

at 

0.6915 

0.0524 

0.8771 

0.6915 

0.0524 

0.8771 

-2.740E-12  -4.007E-11  5.075E-12 

210875 

s  at 

0.2178 

0,5598 

0.0404 

0.2178 

0.5598 

0,0404 

3.556E-10  5,361  E-1 2  2.491  E-11 

211010 

s  at 

0.3193 

0.8268 

0.5607 

0.3193 

0.8268 

0.5607 

1.124E-11  -6.907E-13  •7.970E-12 

211372 

s  at 

0.5618 

0.0342 

0.8184 

0.5618 

0.0342 

0.8184 

-7.687E-12  2.765E-12  2.117E-12 

211560 

s  at 

0.7413 

0.5930 

0.8448 

0.7413 

0.5930 

0.8448 

5.167E-11  1.461E-12  5.340E-13 

211583  X  at 

0.2750 

0.6497 

0.5231 

0.2750 

0.6497 

0.5231 

4.892E-11  1.147E-11  -2.458E-11 

212268 

at 

0.9348 

0.0056 

0.7017 

0.9348 

0.0056 

0.7017 

4.818E-14  3.245E-12  -1.899E-12 

212736 

at 

0.0288 

0.2236 

0.5193 

0.0288 

0.2236 

0.5193 

1,245E-10  6.876E-12  -2.751E-11 

212750 

at 

0,7852 

0.2238 

0.9559 

0,7852 

0.2238 

0.9559 

9.442E-12  6.809E-12  -1.490E-12 

212768 

s  at 

0.5540 

0.2664 

0.1439 

0.5540 

0.2664 

0.1439 

'  -9.761E-12  -8.930E-11  -1.374E-11 

213069 

at 

0.0034 

0.6334 

0.2542 

0.0034 

0.6334 

0.2542 

1.071  E-11  2.340E-13  9.899E-11 

213172 

at 

0.6264 

0.4470 

0.6480 

0,6264 

0.4470 

0.6480 

-2.328E-11  -2.598E-12  -1.183E-11 

213906 

at 

0.0081 

0.3101 

0.0845 

0.0081 

0.3101 

0.0845 

4.887E-11  -1.603E-11  2.651  E-1 2 

214219 

X  at 

0.6005 

0.7319 

0.2025 

0.6005 

0.7319 

0.2025 

-2.264E-12  2.833E-13  -8.709E-12 

214567 

s  at 

0.2668 

0.4752 

0.7056 

0,2668 

0.4752 

0.7056 

6.447E-11  1.114E-10  -1.651E-12 

215047 

at 

0.2101 

0,0830 

0.2754 

0.2101 

0.0830 

0.2754 

-6.411E-12  -1.922E-11  4.824E-11 

216858 

X  at 

0.8672 

0.3129 

0.1632 

0.8672 

0.3129 

0.1632 

1.384E-13  -1.432E-11  -5.370E-12 

217022 

s  at 

0.1488 

0.7041 

0.1598 

0.1488 

0.7041 

0.1598 

-1.077E-11  1.116E-12  -6.309E-12 

217147 

s  at 

0.5464 

0.2553 

0.5463 

0,5464 

0.2553 

0.5463 

-1.230E-11  -1.390E-10  -1.233E-11 

217374 

X  at 

'0.7659 

0.3472 

0.7137 

0.7659 

0.3472 

0.7137 

2.060E-11  -3.729E-12  -1.221E-12 

217636 

at 

0,8302 

0,8488 

0.2350 

0.8302 

0.8488 

0.2350 

1.172E-12  -1.490E-13  1,929E-10 
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219003 

s 

at 

0,1396 

0.3449 

0.1838 

0.1396 

0.3449 

0.1838 

-4.083E-12  -1.881  E-11 

219311 

at 

0.4003 

0.7822 

0,2177 

0,4003 

0.7822 

0.2177 

1.166E-11 

-8.987E-12  3.575E-10 

219529 

_at 

0.0006 

0.1329 

0.2117 

0.0006 

0.1329 

0.2117 

1,365E-12 

4.946E-11  -6.016E-12 

219911 

s 

at 

0.0046 

0.9174 

0.3660 

0.0046 

0.9174 

0.3660 

3.382E-12 

-6.561  E-14  3.267E-11 

219922 

s 

at 

0.6813 

0.8565 

0.2853 

0.6813 

0.8565 

0.2853 

-3.912E-12 

-7.661E-15  3,464E-11 

219948 

X 

at 

0.3751 

0.5325 

0.8177 

0.3751 

0.5325 

0.8177 

2.483E-11 

1.484E-11  2.190E-12 

219978 

s 

at 

0.0000 

0.5662 

0.0158 

0.0000 

0.5662 

0.0158 

1.278E-12 

4.178E-12  3.377E-12 

220330 

s 

at 

0.9907 

0.5827 

0.6795 

0.9907 

0.5827 

0.6795 

2.195E-12  -4.165E-12 

220418 

at 

0.7624 

0.6845 

0,4561 

0.7624 

0.6845 

0.4561 

2.361  E-11 

2.687E-12  -1.658E-10 

220785 

at 

0.4850 

0.1540 

0.1907 

0.4850 

0.1540 

0.1907 

-7.385E-11 

1.468E-11  -1.408E-11 

221586 

s 

at 

0.0726 

0,2021 

0.1936 

0.0726 

0.2021 

0.1936 

6.538E-12 

1.857E-11  ■1.248E-11 

222299 

X 

at 

0.0269 

0.8243 

0.6629 

0.0269 

0.8243 

0.6629 

1.695E-10 

-8.288E-13  -7.283E-12 

Table  5B.  Summary  of  Data 


Program 

Factor 

Lowest 

P  values 

Highest  Mean 

Standard 

deviation 

JAVA 

COND 

2.60E-06 

0.9907 

0.4222 

0.2973 

TIME 

5.58E-03 

0.9718 

0.3873 

0.2839 

TIME*COND 

8.82E-05 

0.9837 

0.4353 

0.2991 

SPSS 

COND 

2.60E-06 

0.9907 

0,4222 

0,2973 

TIME 

5.58E-03 

0.9718 

0.3873 

0.2839 

TIME'COND 

8.82E-06 

0.9837 

0.4353 

0.2991 

Difference 

COND 

-1.79E-10 

3.56E-10 

1.08E-11 

6.40E-11 

TIME 

-2.72E-10 

2.13E-10 

-1.90E-12 

6.08E-11 

TIME*COND 

-1 .66E-10 

3.57E-10 

2.85E-12 

5.50E-11 

COMPARISON  OF  THE  OUTPUT  GENERATED  BY  THE  JAVA  PROGRAM  TO 
OUTPUT  GENERATED  BY  THE  “R”  PROGRAMMING  LANGUAGE  2.1.1 

To  further  validate  the  output  of  the  ANOVA  Java  program,  we  compared  the 
output  it  generated  on  a  full  dataset  (22,283  sequences)  to  that  generated  by  a  routine 
written  in  the  R  programming  language  (v  2.1.1).  As  before,  the  computations  were  run 
on  different  computers,  with  1 5  significant  digits  of  output  generated  (1 6  digits  total).  As 
shown  in  Table  6,  the  differences  between  the  outputs  of  the  two  programs  were  again 
negligibly  small  and  attributable  to  floating  point  errors. 

Importantly,  the  processing  time  required  by  the  two  programs  running  on  the 
same  computer  was  markedly  different.  On  a  late-model  computer  running  a  2.8  GHz 
Pentium  D  dual-core  processor  under  Windows  XP,  the  ANOVA  Java  program  required 
less  than  60  seconds  to  compute  the  results  of  the  dataset  (which  contained  a  total  of 
802,188  data  points).  By  contrast,  the  R  routine,  which  was  designed  only  to  compute 
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the  P  values  for  2-way,  RM  ANOVA  (and  not  1-way  ANOVA  on  time  slopes  or  Tukey’s 
HSD),  required  about  6  %  hours  to  execute. 

On  rare  occasions,  the  ANOVA  Java  program  generated  P  values  of  less  than 
zero.  Of  22,283  sequences  analyzed  in  the  test  dataset,  six  had  P  values  for  COND  of 
less  than  zero,  the  smallest  value  being  -3.09  x  10'^.  In  the  R  programming  language, 
the  corresponding  P  values  for  these  six  sequences  were  all  positive  but  smaller  than  5 
X  1 0'®.  P  values  this  small  are,  for  practical  purposes,  highly  statistically  significant  and 
no  different  from  zero. 


Table  6.  Summary  of  a  comparison  between  the  outputs  of  the  JAVA  ANOVA  program 

to  the  results  generated  by  R  2.1 .1 


Program 

Factor 

Lowest 

P  ’ 

Highest 

tralues 

Mean 

Standard 

deviation 

JAVA 

COND 

-3.08530E-07 

0.99992 

0.37032 

0.31364 

TIME 

1.97614E-08 

1.00000 

0.41854 

0.30518 

COND’TIME 

2.64314E-09 

0.99995 

0.42541 

0.30594 

R 

COND 

3.54267E-12 

0.99992 

0.37032 

0.31364 

TIME 

1.96174E-08 

1.00000 

0.41854 

0.30518 

COND*TIME 

4.44393E-09 

0.99995 

0.42541 

0.30594 

Difference 

COND 

-3.08668E-07 

7.0731 1E-06 

3.74096E-10 

4.82288E-08 

TIME 

-3.31890E-10 

2.81597E-10 

6.27451  E-1 3 

5.43576E-11 

COND'TIME 

-1.80080E-09 

3.23315E-10 

-4.66212E-13 

6.80459E-11 
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DISCUSSION 


The  ANOVA  Java  tool  described  herein  was  developed  because  of  limitations 
inherent  to  popular  general  statistical  packages  that  made  them  unattractive  choices  for 
performing  two-way  ANOVA  on  Affymetrix  microarray  data.  These  limitations  included 
cost  (SAS),  speed  of  execution  (R  programming  language),  and  user  labor  (SPSS).  By 
contrast,  the  ANOVA  Java  tool  allows  the  user  to  rapidly  analyze  any  number  of 
instances,  limited  only  by  the  amount  of  physical  and  virtual  memory  available  on  the 
user's  computer.  Furthermore,  the  ANOVA  Java  tool  is  very  fast  -  it  was  able  to  perform 
a  computation  on  a  dataset  at  a  speed  that  was  two  orders  of  magnitude  faster  than  a 
routine  written  in  R.  The  output  files  generated  by  the  tool  are  tab-delimited  files  that 
can  readily  be  imported  into  almost  any  major  relational  database  management  system, 
such  as  Oracle,  SQL  server,  Microsoft  Access,  etc.,  thus  facilitating  output 
management.  Finally,  the  tool  gives  the  user  the  option  to  determine  the  precision  to 
which  each  P  value  is  computed,  though  for  practical  purposes,  most  will  find  five  to  six 
digits  sufficient. 

At  present,  the  ANOVA  Java  tool  only  allows  one-way  ANOVA  on  slopes  and 
two-way  repeated  measures  ANOVA.  However,  the  data  structure  it  uses  is  very 
generic  and  should  be  readily  capable  of  being  expanded  to  allow  performance  of  multi¬ 
factorial  ANOVA.  Furthermore,  since  the  ANOVA  tool  is  entirely  implemented  in  Java, 
it  can  take  advantage  of  the  platform-independent  database  connection  mechanism 
provided  by  Java  through  JDBC  (7).  Using  JDBC  the  tool  can  be  expanded  to  read  input 
data  from  the  database  and  store  the  results  directly  to  the  database. 

An  additional  limitation  to  the  ANOVA  Java  tool  is  that  it  does  not  perform 
assumption  checking  on  the  data.  This  is  more  of  a  theoretical  than  practical  limitation 
at  present.  ANOVA  methods  in  general  (and  paired,  repeated-measures  ANOVA 
methods  in  particular)  are  fairly  robust  with  regard  to  deviations  from  underlying 
assumptions  methods.  Additionally,  when  the  numbers  of  replicates  are  low,  the  ability 
of  existing  statistical  methods  to  reliably  detect  deviations  from  normality  or 
homogeneity  of  variance  is  often  quite  low. 

Another  important  limitation  to  the  tool  is  that  it  was  designed  with  a  specific 
experiment  in  mind,  and  some  features  of  the  program  (such  as  the  time  interval 
denominator  used  to  compute  slopes  in  the  one-way  ANOVA  routine)  can  only  be 
modified  by  direct  modification  of  the  source  code. 

In  summary,  we  have  designed  and  implemented  a  platform-independent 
software  tool  that  can  rapidly  and  accurately  compute  the  results  of  multifactorial, 
repeated-measures  microarray  experiments.  The  program  accepts  text  input  and 
generates  output  that  can  easily  be  imported  into  popular  commercial  databases,  with  a 
degree  of  precision  that  is  user-defined.  The  tool  thus  overcomes  several  important 
limitations  that  are  evident  when  one  attempts  to  use  general-purpose  software 
packages  to  perform  ANOVA  on  microarray  data. 
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APPENDIX  A 


ILLUSTRATIVE  EXAMPLE  OF  CONFIGURATION  FILE  FORMAT 

[config] 

#  Number  of  factors  to  consider  in  ANOVA. 

NumOfFactors  =  3 

#  Character  that  separates  the  columns  in  the  input  file. 

#  The  value  of  this  parameter  must  be  enclosed  in  “  and 
ColumnSeparator  = 

#  Input  file. 

DataFile  =  input.txt 

TwoWayOutputFile  =  SampleTwoWayAnova.txt 
OneWayOutputFile  =  SampleOneWayAnova.txt 
HeadersForOneWayAnova  =  Factor2 

#  Now  describe  each  factor  individually. 

#  A  factor  paragraph  has  the  format  FactorW,  where  N  is  substituted  by  an 

#  integer. 

#  The  number  of  Factor  paragraphs  must  be  equal  to  the  value  of  NumOfFactors 

#  parameter  in  the  config  paragraph  (see  above). 

#  Factor  "sample"  has  6  experiments,  and  it  is  specified  in  the  second  column  of 

#  the  input  file. 

[Factorl] 

Name  =  sample 
Columns  =  2 
MaxNumOfIndex  =  6 

#  Example  of  a  factor  that  has  categorical  values. 

#  The  factor  is  enumerated  by  assigning  an  integer  to  each  category. 

#  The  integers  must  start  at  1  and  increment  by  1 . 

#  Note  that  the  category  name  and  the  corresponding  integer  are  separated  by  a 

#  tab  character. 

#  For  example,  “heat"  and  “1"  are  separated  by  a  tab. 

[Factor2] 

Name  =  condition 
Columns  =  3 
MaxNumOfIndex  =  3 

Enumerate  =  heat  1  cold  2  control  3 

[Factors] 

Name  =  time 
Columns  =  4 
MaxNumOfIndex  =  2 

#  Note  that  time  here  is  treated  as  a  categorical  parameter,  with  two  values, 

#  1  and  2,  corresponding  to  T=0  and  T=3,  respectively.  The  time  interval  denominator 

#  used  to  compute  slopes  is  presently  not  user-modifiable 


21 


ILLUSTRATIVE  EXAMPLE  OF  INPUT  FILE  FORMAT 

For  the  sample  dataset,  the  experiment  is  performed  6  times  for  3  conditions.  In 
each  trial,  each  gene  is  measured  at  two  diTerenl  time  points  in  every  condition. 
Therefore,  we  will  have  36  values  for  each  gene.  The  matrix  for  one  gene  is  described 
in  Table  1 .  However,  in  the  matrix  each  row  provides  measured  intensities  of  the  gene 
at  two  different  time  points  for  a  given  sample  and  a  given  condition.  The  matrix  in  Table 
1  is  essentially  a  3-dimensional  matrix  of  size  3  x  6  x  2.  It  is  a  representation  of  a  3- 
dimensional  matrix  in  2-dimension  space.  Note  that  the  matrix  in  Table  1  contains  3 
embedded  smaller  tables,  one  for  each  condition.  The  problem  of  representing  the  data 
in  matrix  form  becomes  difficult  as  we  Increase  the  number  of  factors  {i.e.,  the 
dimensions  of  the  matrix). 

However,  the  matrix  in  Table  1  can  be  represented  in  a  different  format  with  N  + 

1  columns,  where  the  first  N  columns  contain  indices  for  the  number  of  experiment,  the 
condition,  and  the  time  point,  respectively,  and  the  last  column  contains  the  expression 
value  of  the  gene.  This  format  can  describe  one  gene  only.  If  we  add  the  unique  ID  of 
each  gene  as  the  first  column,  we  can  now  represent  the  data  for  any  number  of  genes. 
The  new  table  now  contains  N  +  2  columns.  Table  7  shows  the  format  with  N  +  2 
columns  separated  by  commas,  where  the  'irst  column  is  the  unique  ID  of  the  gene,  the 
second  column  is  the  number  of  the  experiment,  the  third  column  is  the  condition,  the 
fourth  column  is  the  time  (treated  as  a  categorical  variable  with  1  corresponding  to  the 
first  time  point  and  2  corresponding  to  the  second  time  point),  and  the  last  column  is  the 
expression  value  of  the  gene.  Note  that  this  format  requires  that  each  instance  in  the 
input  file  be  assigned  a  unique  label  even  if  the  input  file  contains  only  one  instance. 

Table  7.  Sample  Input  file. 


AFFX-BioB-5_at,'  ,Control,  1 ,470.90 
AFFX-BioB-5_at,', Control, 2, 455.50 
AFFX-BioB-5_at,' , Heat, 1 ,486.70 
AFFX-BioB-5_at,'  ,Heat.2, 464.40 
AFFX-BioB-5_at,'  ,Cold,1 ,378.40 
AFFX-BioB-5_at,-  ,Cold,2,380.00 
AFFX-BioB-5_at,2,Controi,1, 370.40 
AFFX-BioB-5_at,2,Control,2,399.30 
AFFX-BioB-5_at,2,Heat,1 ,41 8.20 
AFFX-BioB-5_at, 2, Heat,2, 454.70 
AFFX-BioB-5_at,2, Cold, 1,395.30 
AFFX-BioB-5_at,2,Cold,2,425.40 
AFFX-BioB-5_at,3, Control, 1,393.40 
AFFX-BioB-5_at,3,Control, 2, 392.70 
AFFX-BioB-5_at,3,Heat,1 ,386.60 
AFFX-BioB-5_at, 3, Heat,2, 428.60 
AFFX-BioB-5_at,3, Cold, 1,282.90 
AFFX-BioB-5_at,3,Cold,2,331 .80 
AFFX-BioB-5  at, Control, 1,347.80 
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AFFX-BioB-5_at,4.Control, 2, 389.80 
AFFX-BioB-5_at,4,Heat,1 ,393.80 
AFFX-BioB-5_at,4,Heat,2,443.60 
AFFX-BioB-5_at,4,Cold,1 ,354.20 
AFFX-BioB-5_at, 4, Cold, 2,432.30 
AFFX-BioB-5_at, 5, Control, 1,339.50 
AFFX-BioB-5_at,5,Control, 2,373.10 
AFFX-BioB-5_at, 5, Heat, 1,402.40 
AFFX-BioB-5_at, 5, Heat,2, 524.00 
AFFX-BioB-5_at,5,Cold,  1 ,346.1 0 
AFFX-BioB-5_at,5,Coid,2, 360.60 
AFFX-BioB-5_at,6,Control, 1,354.50 
AFFX-BioB-5_at,6, Control, 2, 31 4.40 
AFFX-BioB-5_at,6,Heat,1, 376.80 
AFFX-BioB-5_at,6,Heat,2, 1 097. 1 0 
AFFX-BioB-5_at,6,Co]d,  1 ,484.50 
AFFX-BioB-5  at,6,Cold,2,419.2Q 


ILLUSTRATIVE  EXAMPLE  OF  OUTPUT  FILE  FORMAT  FOR  TWO-WAY  ANOVA 
WITH  REPEATED  MEASURES  ON  ONE  FACTOR 

An  example  of  the  output  generated  by  the  two-way  ANOVA  algorithm  Is  listed  in  Table 
8.  The  program  generates  p-values  for  main  effect  (condition),  the  repeated  measures 
factor  (time),  and  the  interaction  of  condition  x  time. 

Table  8.  Sample  output  file  for  two-way  ANOVA  with  repeated  measures  on  one  factor. 


Gene  Name 

p-value  for  condition 

p-value  for  time 

p-value  for  interaction 

AFFX-BioC-5  at 

0,242995 

0.301010 

0.525759 

AFFX-Crex-3  at 

0.135440 

0.083110 

0.429870 

AFFX-PheX-M  at 

0.571269 

0.881806 

0.122794 

AFFX-DapX-M  at 

0.940233 

0.689611 

0.081855 

AFFX-ThrX-M  at 

0.704193 

0.595223 

0.769351 

AFFX-LysX-5  at 

0.844022 

0.777465 

0.868368 
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ILLUSTRATIVE  EXAMPLE  OF  OUTPUT  FILE  FORMAT  FOR  ONE-WAY  ANOVA 

The  sample  output  for  one-way  ANOVA  is  shown  in  Table  9.  The  tool  provides 
confidence  intervals  for  three  different  p-values:  0.01 , 0.05,  and  0.1 .  The  order  in  which 
the  columns  are  presented  is  as  follows: 

Column  1  Sequence  name 

Column  2  P  value  for  condition 

Subsequent  columns:  Tukey  HSD  confidence  intervals  for  subgroup  comparisons  at  P- 
values  0.01 , 0.05,  and  0.1 .  The  number  of  columns  will  depend  upon  the  number 
of  subgroups  being  compared,  as  define  by  “Factor2”  in  the  configuration  file.  In 
our  sample  computation,  there  are  three  subgroups  (heat,  cold,  control)  and 
hence,  three  possible  comparisons  (heat  vs.  cold,  heat  vs.  control,  cold  vs. 
control).  So,  in  our  output  example: 

Columns  3,  4,  5:  Confidence  intervals.  Index  (1)  vs.  Index  (2)  (in  this  case.  Heat 
vs.  Cold)  for  the  following  P  values: 

Column  3:  P=0.01 
Column  4:  P=0.05 
Column  5:  P=0.10 

Columns  6,  7,  8:  Confidence  intervals.  Index  (1)  vs.  Index  (3)  (in  this  case.  Heat 
vs.  Control) 

Column  6:  P=0.01 
Column  7:  P=0.05 
Columns:  P=0.10 

Columns  9,  10,  11:  Confidence  intervals.  Index  (2)  vs.  Index  (3)  (in  this  case. 
Cold  vs.  Control)  for  the  following  P  values: 

Column  9:  P=0.01 
Column  10:  P=0.05 
Column  11:  P=0.10 
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Table  9.  Illustrative  example  of  output  file  format  for  one-way  ANOVA. 
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APPENDIX  B 


PROGRAM  LISTING 

The  java  code  to  run  an  ANOVA  analysis  has  been  laid  out  in  two  packages,  namely, 
“anova. client”  and  “anova.util”.  The  package  “anova. client”  contains  one  file  named 
“AnovaAnalysisI  Java”,  which  is  the  main  application.  The  list  of  the  files  contained  in 


the  package  “anova.util”,  is  as  follows: 

1. 

DataFile.java 

2. 

Factorjava 

3. 

inlElement.java 

4. 

IniFile.java 

5. 

IniParagraph.java 

6. 

InputFile.java 

7. 

OutputFile.java 

The  application  also  uses  an  external  library  named  “commons-math-1 .0-RC1  Jar”, 
available  from  apache  [7].  The  actual  code  contained  in  the  files  listed  above  is  as 
follows: 

AnovaAnalysisI  .lava 
/• 

*  Created  on  Aug  28,  2003 
*1 

package  anova. client; 

r 

*  ©author  Purvesh  Khatri 
*/ 

import  java.io.FileWriter; 
import  Java. io.lOException; 
import  java.io.PrintWriter; 
import  java. util. Date; 
import  java. util. Iterator; 
import  java. util. StringTokenizer; 

import  org.apache.commons.math.MathException; 

import  org. apache. commons.math. distribution. DistributionFactory; 

import  org. apache. commons. math. distribution. FDistribution; 

import  anova.util. DataFile; 
import  anova.util. Factor; 
import  anova.util. IniFile; 

public  class  AnovaAnalysisI 
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{ 

private  static  int  numOfSignificantDigits  =  8; 

public  static  void  main(String  Q  args) 

{ 

if  (args. length  <  1) 

{ 

System. err.printlnC'USAGE:  java  anova  client.AnovaAnalysis  <config_file>" 
+  "  [-d  <number  of  significant  digits>]'  ): 

System. err.printin( 

"Also  make  sure  commons-math  library  is  in  classpath."); 

System. err.println( 

"You  may  need  to  modify  the  classpath  using  -cp  argument."); 

System. exit(O); 

} 

IniFile  ini  =  null; 

//  Read  the  configuration  file. 

try 

{ 

ini  =  new  lniFile(args[0]); 

} 

catch  (lOException  e) 

{ 

System. out.println(e.getMessage()); 

System. exit(1); 

} 

try 

{ 

if  (args.lenglh  >  1  &&  args[1].equals("-d*)) 

{ 

AnovaAnalysisl. numOfSignificantDigits  =  lnteger.parselnt(args[2],  10); 

} 

} 

catch  (NumberFormatException  nfe) 

{ 

nfe.printStackTraceO: 

System. err.println{"NumberFormatException:  Invalid  number  of " 

+  "significant  digits  specified: "  +  args;2]); 

System. err.println("\tlt  must  be  an  integer  number."); 

System. err.println("Using  default  number  of  significant  digits."); 
AnovaAnalysislnumOfSignificantDigits  =  8; 

} 
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//  Create  an  object  to  parse  the  input  file,  based  on  configuration  file. 
DataFile  dataFile  =  new  DataFile{ini): 

//  Now  actually  parse  the  file. 

Date  start  =  new  Date(): 

System. err.print("Reading  file..."); 
dataFile. readDataFileQ: 

Date  readFile  =  new  Date{): 

System. err.println("Done.  Time  taken: " 

+  ((readFile.getTimeO  -  start.getTimeO)  / 1000)  +  "  sec"); 

if  (  !ini.getProfile("config",  "twowayoutputfile").equals("")) 
doTwoWayAnalysis{dataFile,  ini,  dataFile.getSizeArrayO); 

Date  endTwoWay  =  new  Date(); 

if  ( !ini.getProfile{"config",  "twowayoutputfile").equals("")) 

System. err.printlnC’Two  way  ANOVA  done.  Time  for  analysis: " 

+  ({endTwoWay.getTimeO  -  readFile.getTimeO)  / 1000)  +  "  sec"); 

if  {  lini.getProfileC'config",  "onewayoutputfile").equals("")) 
doOneWayAnalysis(dataFile,  ini,  dataFile.getSizeArrayO,  true); 

Date  endOneWay  =  new  DateO; 

if  (  lini.getProfileC'config",  "onewayoutputfile").equats("")) 

System. err.printinC'One  way  ANOVA  done.  Time  for  analysis: " 

+  ((endOneWay.getTimeO  -  endTwoWay.getTimeO)  / 1000)  + "  sec"); 

System. err.printlnC’  Total  time: " 

+  {(endOneWay.getTimeO  -  start.getTimeO)  / 1000)  +  "  sec"); 

} 

I** 

*  @param  value 

*  Floating-point  number  that  needs  to  be  rounded  off 

*  @param  places 

*  Number  of  significant  digits  to  which  value  should  be  rounded  off 

*  ©return  Returns  the  rounded  off  number. 

*/ 

public  static  double  round(double  value,  int  places) 

{ 

long  factor  =  (long)  Math.pow(10,  places): 

//  Shift  the  decimal  the  correct  number  of  places 
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//  to  the  right, 
value  =  value  *  factor; 

//  Round  to  the  nearest  integer, 
long  tmp  =  Math.round(value); 

//  Shift  the  decimal  the  correct  number  of  places 

//  back  to  the  left. 

return  (double)  tmp  /  factor; 

} 

private  static  void  printHeadersfPrintWriter  pw,  IniFile  iniFile) 

{ 

String  headerFactor  =  iniFile. getProfile("config",  "HeadersForOneWayAnova"); 
String  headersLine  =  iniFile. getProfile(headerFactor,  "enumerate"); 
if  (headersLine  ==  null  ||  headersLine. equals("")) 
throw  new  Nu!IPointerException(headerFactor  + "  must  define  headers " 

+  "using  Enumerate"); 

StringTokenizer  stk  =  new  StringTokenizer(headersLine.substring( 
headersLine. indexOf(-')  +  1)); 

//  the  headersLine  must  be  tokenized  in  even  number  of  tokens 
String  Q  headers  =  new  String  [stk.countTokensO  /  2]: 
for  (int  i  =  0;  i  <  headers. length;  i++) 

{ 

headers[i]  =  stk.nextToken().trim(): 
stk.nextTokenO; 

} 

for  (int  i  =  0;  i  <  headers. length;  i++) 

{ 

for  (int  j  =  i  +  1 ;  j  <  headers. length;  j++) 

{ 

II  The  order  of  p-values  is  hard-codes  because  we  will  always  calculate 
II  confidence  intervals  in  this  order. 

//  see  doTukeyComparison(). 

pw.print('T  +  headers[i]  +  "  - "  +  headersO]  +  "  (p  =  0.01 )"); 

pw.print("\t"  +  headers[i]  +  "  - "  +  headersO]  "  (P  =  0.05)"); 

pw.print("\t"  +  headersp]  +  "  - "  +  headersO]  "  (p  =  0.1)"); 

} 

} 

} 

private  static  void  doOneWayAnalysis(DataFile  dataFile,  IniFile  iniFile, 
int  []  sizeArray,  boolean  doTukeyComparisons) 


30 


{ 

PrintWriter  pw  =  null; 
try 
{ 

pw  =  new  PrintWriter(new  FileWriter( 
iniFile.getProfile{"config",  "onewayoutputfile''))); 
pw.print("Name\tp  value  for  condition"); 

if  (doTukeyComparisons) 
printHeaders(pw,  iniFile); 

pw.printlnO; 

} 

catch  (lOException  ioe) 

{ 

System. err.println(ioe): 

System. exit(4): 

} 

Iterator  instancelterator  =  dataFile.getKeySetIteratorO: 

int  countConditionPValue  =  0; 

while  (instancelterator. hasNextO) 

{ 

II  Get  the  next  instance. 

String  name  =  (String)  instancelterator.next(): 

Factor  factor  =  dataFile.get(name); 

int  Q  indexArray  =  new  int  [3]; 

int  i  =  0,  j  =  0; 

II  calculate  means  for  each  conditions. 

double  Q  conditionMean  =  new  double  [si2eArray[1]]:  II  Yi. 

for  (j  =  0;  j  <  sizeArray[1];  j++) 

{ 

for  (i  =  0;  i  <  sizeArray[0]:  i++) 

{ 

indexArray[0]  =  i; 
indexArray[1]  =  j; 
indexArray[2]  =  0; 

double  f1  =  ((Double)  factor.get(indexArray)).doubleValue(): 
indexArray[2]  =  1 ; 

double  f2  =  ((Double)  factor.get(indexArray)).doubleValue(): 
conditionMeanO]  +=  ((f2  -  f1 )  /  (3)); 


31 


} 

conditionMeanO]  =  conditionMeanO]  /  I'i): 

} 

//  calculate  global  mean 
double  globalMean  =  0.0;  //  Y.. 
for  0  =  0;  j  <  sizeArray[1]:  j++) 

{ 

for  (i  =  0;  i  <  sizeArray[0]:  i++) 

{ 

indexArray[0]  =  i; 
indexArray[1]  =  j; 
indexArray[2]  =  0; 

double  f1  =  ((Double)  factor.get(indexArray)).doubleValue(): 
indexArray[2]  =  1; 

double  f2  =  ((Double)  factor.get(indeKArray)).doubleValue(): 
globalMean  +=  ((f2  -  f1 )  /  (3)); 

} 

} 

globalMean  =  globalMean  /  (I  *  j); 

double  mstr  =  0.0; 

for  (j  =  0;  j  <  siz6Array[1]:  j++) 

{ 

mstr  +=  (Math.pow(conditionMean[j]  -  globalMean,  2)  *  sizeArray[0]); 

} 

mstr  =  mstr  /  (sizeArray[1]  - 1 ): 
double  mse  =  0.0; 

for  (j  =  0;  j  <  sizeArray[1]:  j++)  //  For  each  condition 

{ 

for  (i  =  0;  i  <  slzeArray[0];  i++)  //  each  measurement  in  each  condition 

{ 

indexArray[0]  =  I; 
indexArray[1]  =  j: 
indexArray[2]  =  0; 

double  f1  =  ((Double)  factor.get(inde)(Array)).doubieValue(); 
indexArray[2]  =  1; 

double  f2  =  ((Double)  factor.get(inde.<Array)).doubleValue(); 
mse  +=  Math.pow(((f2  -  f1 )  /  (3))  -  conditionMean[j],  2); 

} 

} 

mse  =  mse  /  ((sizeArray[0]  *  sizeArray[1])  -  sizeArray[1]): 


try 

{ 
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//  Now  calculate  p-values  for  condition. 

DistributionFactory  factory  =  DistributionFactory.newlnstance(); 

FDistribution  fDist  =  factory.createFDistribution(sizeArray[1]  - 1 , 
(sizeArray[0]  *  sizeArray[1])  -  sizeArray[1]); 
double  fStar  =  mstr  /  mse; 

double  conditionPValue  =  1  -  fDist.cumulativeProbability{fStar); 
if  (conditionPValue  <=  0.01 ) 
countConditionPValue++: 

pw.print(name  + 

if  (conditionPValue  <=  1e-10) 
pw.print(O.O): 
else 

pw.print(round(conditionPValue, 

AnovaAnalysisl.numOfSignificantDigits)): 

} 

catch  (MathException  me) 

{ 

System. err.prinlln("MathException  for  instance  "  +  name); 
me.printStackTraceO; 

System.err.prinlln("Skipping  the  instance "  +  name  +  "..."); 

} 

if  (doTukeyComparisons) 

{ 

doTukeyComparison(conditionMean,  mse,  sizeArray,  pw); 

} 

pw.printInO: 

}  //  end  of  while  loop. 

pw.closeO: 

} 

private  static  void  doTukeyComparison(double  Q  conditionMean,  double  mse, 
int  □  sizeArray,  PrintWriter  pw) 

{ 

double  D  g  =  {4.83,  3.67,  3.14}; 
int  i,j,  k; 

for  (i  =  0;  i  <  conditionMean. length;  i++) 

{ 

for  (j  =  i  +  1;  j  <  conditionMean. length;  j++) 

{ 
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for  (k  =  0;  k  <  g. length;  k++) 

{ 

double  approxD  =  conditionMean[i]  -  oonditionMean]]]; 
double  approxStdD  =  Math.sqrt( 

mse  *  (1  /  (double)  sizeArray[0]  +  1  /  (double)  sizeArray[0])): 
double  t  =  g[k]  /  Math.sqrt(2); 

double  lowerBound  =  approxD  -  (t  *  approxStdD); 
double  upperBound  =  approxD  +  (t  *  approxStdD); 

pw.print("\t['' 

+  round(lowerBound,  AnovaAnalysisI  .numOfSignificantDigits)  +  ", " 
+  round(upperBound,  AnovaAnalysisI. numOfSignificantDigits)  +  "]"); 

} 

} 

} 

} 

private  static  void  doTwoWayAnalysis(DataFiie  dataFile,  IniFile  iniFile, 
int  D  sizeArray) 

{ 

PrintWriter  pw  =  null; 

try 

{ 

pw  =  new  PrintWriter(new  FileWriter( 

iniFile. getProfile("config",  "twowayoutputfile"))); 
pw.printtn("Probe\tp  value  for  conditionMp  value  for  time\t " 

+  "p  value  for  interaction"); 

} 

catch  (lOException  ioe) 

{ 

System. err.println(ioe): 

System. exit(4); 

} 

Iterator  instancelterator  =  dataFile. getKeySetlterator(); 

int  countConditionPValue  =  0; 
int  countTimePVaiue  =  0; 
int  countInteractionPValue  =  0; 

while  (instancelterator. hasNextO) 

{ 

String  name  =  (String)  instancelterator.rext(); 

Factor  factor  =  dataFile. get(name); 
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int  n  indexArray  =  new  int  [3]; 
int  i  =  0,  j  =  0,  k  =  0; 


II  calculate  means  for  each  conditions. 

double  Q  conditionMean  =  new  double  [slzeArray[1]];  II  Y.j. 

for  (j  =  0;  j  <  sizeArray[1];  j++) 

{ 

for  (i  =  0;  i  <  sizeArray[0]:  I++) 

{ 

for  (k  =  0;  k  <  sizeArray[21:  k++) 

{ 

indexArray[0]  =  i; 
indexArray[1]  =  j; 
indexArray[2]  =  k; 

conditionMean[j]  +=  ((Double)  factor.get(indexArray)).doubleValue(}; 

} 

} 

conditionMeanOl  =  conditionMean|j]  /  (k  *  I); 

} 

//  calculate  global  mean 
double  globalMean  =  0.0;  II Y... 
for  G  =  0:  j  <  sizeArray[1];  j++) 

{ 

for  (I  =  0;  i  <  sizeArray[0]:  i++) 

{ 

for  (k  =  0;  k  <  sizeArray[2];  k++) 

{ 

indexArray[0]  =  i; 
indexArray[1]  =  j; 
indexArray[2]  =  k; 

globalMean  +=  ((Double)  factor.get(lndexArray)).doubleValue(); 

} 

} 

} 

globalMean  =  globalMean  /  (i  *  j  *  k); 

II  calculate  means  for  each  time  point. 

double  Q  timeMean  =  new  double  [slzeArray[2]];  //  Y..k 

for  (k  =  0;  k  <  slzeArray[2]:  k++) 

{ 

for  (i  =  0;  i  <  sizeArray[0]:  i++) 

{ 

for  (j  =  0;  j  <  sizeArray[1]:  j++) 

{ 
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indexArray[0]  =  i; 
indexArray[1]  =  j; 
indexArray[2]  =  k; 

timeMean[k]  +=  ((Double)  factor.get(indexArray)).doubleValue(); 

} 

} 

timeMean[k]  =  timeMean[k]  /  (j  *  0: 

} 

//  calculate  contime  time  interaction  mean. 

//  Y.jk 

double  QQ  conditionTimeMean  =  new  double[sizeArray[1]][sizeArray[2]]: 
for  (k  =  0;  k  <  sizeArray[2];  k++) 

{ 

for  (j  =  0;  j  <  sizeArray[1]:  j++) 

{ 

for  (i  =  0;  i  <  sizeArray[0]:  i++) 

{ 

indexArray[0]  =  i; 
indexArray[1]  =  j; 
indexArray[2]  =  k; 
conditionTimeMean[j][k]  += 

((Double)  factor.get(indexArray)).doubleValue(): 

} 

conditionTimeMean[j][k]  =  conditionTimeMeanOKk]  /  (i); 

} 

} 

//Yij. 

double  QD  sampleConditionMean  =  new  double  [sizeArrayfO]]  [sizeArray[1]]: 
for  (i  =  0;  i  <  slzeArray[0]:  i++) 

{ 

for  (j  =  0;  j  <  slzeArray[1];  j++) 

{ 

for  (k  =  0;  k  <  sizeArray[2]:  k++) 

{ 

indexArray[0]  =  I; 
indexArray[1]  =  j; 
indexArray[2]  =  k; 
sampleConditionMean[i]0]  += 

((Double)  factor.get(lndexArray)).doubleValue(): 

} 

sampleConditionMean[i]0]  =  sampleConditionMean[i]0]  /  (k); 

} 

} 
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//  Now  calculate  mean  sum  of  squares. 

double  msa  =  0.0; 

for  (j  =  0;  j  <  sizeArray[1];  j++) 

{ 

msa  +=  Math.pow(conditionMeanO]  -  globalMean,  2); 

} 

msa  =  (msa  *  si2eArray[2]  *  sizeArray[0])  /  (sizeArray[1  ]  - 1); 

double  msb  =  0.0; 

for  (k  =  0;  k  <  sizeArray[2];  k++) 

{ 

msb  +=  Math.pow(timeMean[k]  -  globalMean,  2); 

} 

msb  =  (msb  *  sizeArray[1]  *  sizeArray[0])  /  (sizeArray[2]  - 1); 

double  msab  =  0.0; 

for  (j  =  0:  j  <  sizeArray[1];  j++) 

{ 

for  (k  =  0;  k  <  sizeArray[2];  k++) 

{ 

msab  +=  Math.pow(conditionTimeMean[j][k]  -  conditionMeanO] 

-  timeMean[k]  +  globalMean,  2); 

} 

} 

msab  =  (msab  *  sizeArray[0])  /  ((sizeArray[1]  - 1)  *  (sizeAnray[2]  - 1)); 

double  mssa  =  0.0; 

for  (I  =  0;  i  <  sizeArray[0]:  i++) 

{ 

for  (j  =  0;  j  <  sizeArray[1];  j++) 

{ 

mssa  +=  Math.pow(sampleConditionMean[i]0]  -  conditionMean[j],  2); 

} 

} 

mssa  =  (mssa  *  sizeArray[2])/  (sizeArray[1]  *  (sizeArray[0]  - 1)); 

double  msbsa  =  0.0; 

for  (i  =  0;  i  <  sizeArray[0]:  i++) 

{ 

for  (j  =  0;  j  <  sizeArray[1];  j++) 

{ 

for  (k  =  0;  k  <  sizeArray[2]:  k++) 

{ 

indexArray[0]  =  i; 
indexArray[1]  =  j; 
indexArray[2j  =  k; 
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msbsa  +=  Math.pow(((Double)  factor.get(indexArray)).doubleValue() 
-  conditionTimeMean[j][k]  -  samp,eConditionMean[i]0] 

+  conditionMeanO],  2); 

} 

} 

} 

msbsa  =  msbsa  /  (sizeArray[1]  *  (sizeArray[0]  - 1)  *  (sizeArray[2]  - 1)): 

//  Now  calculate  p-values  for  condition. 

DistributionFactory  factory  =  DistributionFactory.newlnstance(): 

FDistribution  fDist  =  factory.createFDistribution(sizeArray[1]  - 1 , 
sizeArray[1]  *  (sizeArray[0]  - 1 )); 
double  fStar  =  msa  /  mssa: 
double  conditionPValue  =  -1 .0; 
try 
{ 

conditionPValue  =  1  -  fDist.cumulativeProbability(fStar); 
if  (conditionPValue  <=  1  e-1 0) 
conditionPValue  =  0.0; 
if  (conditionPValue  <=  0.01) 
countConditionPValue++; 

} 

catch  (MathException  me) 

{ 

System. err.printlnC'MathException  for  instance  "  +  name); 
me.printStackT  race(); 

System. err.println("Setting  condition  p-value  for  instance  "  +  name 
+  "  as  -1 .0  in  the  output  file."); 
conditionPValue  =  -1; 

} 

fDist.setNumeratorDegreesOfFreedom(sizeArray[2]  - 1); 
fDist.setDenominatorDegreesOfFreedom(sizeArray[1]  *  (sizeArray[0]  - 1) 
*  (sizeArray[2]  - 1 )); 
fStar  =  msb  /  msbsa; 
double  timePValue  =  -1 .0; 
try 
{ 

timePValue  =  1  -fDist.cumulativeProbability(fStar): 
if  (timePValue  <=  1  e-1 0) 
timePValue  =  0.0; 
if  (timePValue  <=  0.01) 
countTimePValue++; 

} 
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catch  (MathException  me) 

{ 

System. err.println("MathException  for  instance  "  +  name): 
me.printStackT  race(): 

System. err.printlnC'Setting  time  p-value  for  instance "  +  name 
+  "  as  -1 .0  in  the  output  file."); 
timePVaiue  =  -1 ; 

} 

fDist.setNumeratorDegreesOfFreedom{{sizeArray[1]  - 1) 

*  (sizeArray[2]  - 1)); 

fDist.setDenominatorDegreesOfFreedom(sizeArray[1]  *  (sizeArray[0]  - 1) 

*  (sizeArray[2]  - 1)); 
fStar  =  msab  /  msbsa; 

double  interactionPValue  =  -1.0; 
try 
{ 

interactionPValue  =  1  -  fDist.cumulativeProbability(fStar): 
if  (interactionPValue  <=  1e-10) 
interactionPValue  =  0.0; 
if  (interactionPValue  <=  0.01) 
countlnteractionPValue++: 

} 

catch  (MathException  me) 

{ 

System. err.printlnC'MathException  for  instance "  +  name); 
me.printStackTraceO; 

System. err.printlnC'Setting  interaction  p-value  for  instance  "  +  name 
+  "  as  -1 .0  in  the  output  file."); 
interactionPValue  =  -1 ; 

} 

pw.println(name  +  "\t" 

+  round(conditionPValue,  AnovaAnalysisI  .numOfSignificantDigits) 

+  "\t"  +  round(timePValue,  AnovaAnalysisI. numOfSignificantDigits) 

+  "\t" 

+  round(interactionPValue,  AnovaAnalysisI  .numOfSignificantDigits)); 


pw.printInC'No.  of  Genes  (P  <=  0.01  )\t"  +  countConditionPValue  +  "\t" 
+  countTimePValue  +  "\t"  +  countlnteractionPValue); 
pw.closeO; 

} 

} 
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DataFile.iava 


r 

*  Created  on  Aug  28,  2003 

*/ 

package  anova.util; 

r 

*  ©author  Purvesh  Khatri 
*1 

import  java.io.*: 
import  java. util.*: 

public  class  DataFile 

{ 

private  IniFile  iniFile; 

private  int  Q  sizeArray; 

private  LinkedHashMap  dataTable; 

public  DataFile(lniFile  ini) 

{ 

iniFile  =  ini; 

dataTable  =  new  LinkedHashMap(); 

} 

r 

*  Reads  the  input  file.  This  is  the  public  interface  provided  by  DataFile 

*  class. 

*/ 

public  void  readDataFile() 

{ 

int  numOfFactors  =  0; 

try 

{ 

numOfFactors  =  lnteger.parselnt(iniFiie.getProfile("config", 
"numoffactors")); 

} 

catch  (NumberFormatException  nfe) 

{ 

System. err.println( 

"DataFile:  Invalid  number  of  factors  specified." 

+  "\nNumOfFactors  must  be  an  integer  number. " 

+  "Please  correct  your  configuration  file."); 

System. err.println("Specified  NumOfFactors: " 
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+  iniFile.getProfileC'config",  "numoffactors")); 

System. exit(1): 

} 

II  find  out  how  the  columns  in  the  file  are  separated. 

String  columnSeparator  =  iniFile.getProfileC’config",  "columnseparator"); 
columnSeparator  = 

columnSeparator.substring(1,  columnSeparatorJengthQ  - 1): 

II  Find  out  the  dimensions  of  the  data  matrix. 
sizeArray  =  createSizeArray(numOfFactors): 

BufferedReader  br  =  null; 

try 

{ 

br  =  new  BufferedReader(new  FileReader( 
iniFile.getProfileC'config",  "datafile"))); 

} 

catch  (lOException  ioe) 

System. err.printlnC’DataFile:  Error  opening  the  datafile."); 

System. err.printlnfPlease  check  file  name: " 

+  iniFile.getProfileC'config",  "datafile")); 

} 

readFile(br); 

} 

private  int  Q  createSizeArray(int  numOfFactors) 

^  List  paragraphNames  =  new  ArrayList(iniFile.getParagraphNames()); 
Collections.sort(paragraphNames); 
int  []  sizeArray  =  new  int  [numOfFactors]; 
intj  =  0; 

int  maxNumOflndex  =  0; 

for  (int  i  =  0;  i  <  paragraphNames.size{);  i++) 

{ 

String  name  =  (String)  paragraphNames.get(i); 
if  (name.toLowerCase().startsWith("factor")) 

{ 

try 

{ 
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maxNumOfIndex  =  lnteger.parselnt( 

iniFile.getProfile(name,  "MaxNumOfIndex")}; 

} 

catch  (NumberFormatException  nfe) 

{ 

System. err.print!n( 

"DataFiie;  Invalid  number  of  indices  specified." 

+  "\nMaxNumOflndex  must  be  an  integer  number. " 

+  "Please  correct  your  configuration  file  " 

+  "for  the  following  paragraph."): 

System. err.println{"Erroneous  Paragraph: "  +  name); 

System. exit(1 ); 

} 

sizeArraylJ]  =  maxNumOfIndex; 
j++: 

} 

} 

return  sizeArray; 

} 

private  void  readFile(BufferedReader  br) 

{ 

String  s  =  null; 

try 

{ 

while  ({s  =  br.readLineO)  !=  null) 

{ 

processlnstance(s): 

} 

} 

catch  (lOException  ioe) 

{ 

System. err.println{"DataFile:  Error  reading  data  file."); 

System. err.println("Error  encountered  when  reading: "  +  s); 

System. exit(2); 

} 

} 

private  void  processlnstance(String  instanceStr) 

{ 

StringTokenizer  stk  =  new  StringTokenizer(instanceStr,  iniFile.getProfile( 
"config",  "ColumnSeparator")); 
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//  first  column  is  always  Instance  id  or  name. 

String  name  =  stk.nextToken(); 

int  columnNum  =  2; 

int  Q  indexArray  =  new  int  [Integer. parse! nt( 
iniFile.getProfileC'config'',  "numoffactors"))]; 

for  (int  i  =  0;  i  <  indexArray.length;  i++) 

{ 

indexArray[i]  =  interpretlndex(stk.nextToken(),  columnNum)  - 1; 
columnNum++; 

} 

Factor  factor  =  (Factor)  dataTable.g6t{name): 

if  (factor  ==  null) 

{ 

factor  =  new  Factor(sizeArray): 

} 

factor.add(indexArray,  new  Doub!e(stk.nextToken().trim())): 
dataTable.put(name,  factor); 

} 

private  int  interpretlndex(String  indexStr,  int  columnNum) 

{ 

List  paragraphNames  =  new  ArrayList(iniFtle.getParagraphNames()): 

Collections.sort(paragraphNames); 

List  factorList  =  new  ArrayList(lnteger.parselnt( 
iniFile.getProfileC'config",  "numoffactors"))); 
for  (int  i  =  0;  i  <  paragraphNames.size();  i++) 

{ 

String  name  =  (String)  paragraphNames.get(i); 
if  (name.toLowerCase().startsWith("factor")) 

{ 

factorList.add(name): 

} 

} 

for  (int  i  =  0;  i  <  factorList.size();  i++) 

{ 

int  column  =  -1; 
try 
{ 
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column  =  lnteger.parselnt(iniFile.getProfile( 

(String)  factorList.get(i),  "Columns"]); 

} 

catch  (NumberFormatException  nfe) 

{ 

System. err.printlnC'DataFile;  Error  parsing  the  column  number  for " 

+  factorList.get(i)): 

System. err.printlnC'Erroneous  line  in  config  file: " 

+  iniFile.getProfile((String)factorLlst.get(i).  "Columns")); 

} 

if  (column  !=  columnNum) 
continue; 
else 
{ 

if  (inlFile.getProfile( 

(String)  factorList.get(i),  "Enumerate"). equals("")) 

{ 

try 

{ 

int  index  =  Integer.parselnt(indexStr); 
return  index; 

} 

catch  (NumberFormatException  nfe) 

{ 

System. err.println( 

"DataFile;  Non-integer  number  specified  in  column  " 

+  columnNum); 

System. exit(3): 

} 

} 

else 

return  interpretEnumeration(indexStr, 

iniFile.getParagraph((String)  factorList.get(i))); 

} 

} 

return  -1; 

} 

private  int  interpretEnumeration(String  indexStr,  iniParagraph  paragraph) 

{ 

String  inputStr  =  indexStr.toLowerCase(]: 

String  enumStr  =  paragraph.tagValue("Enumerate"): 

StringTokenizer  stk  =  new  StringTokenizer(enumStr, "  \t"): 
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ini  enumValue  =  -1; 
while  (stk.hasMoreTokensO) 

{ 

String  enumName  =  stk.nextToken().trim().toLowerCase(); 

try 

{ 

enumValue  =  lnteger.parselnt(stk.nextToken{).trim()): 

} 

catch  (NumberFormatException  nfe) 

{ 

System. err.println("DataFile:  Error  interpreting  enumeration"); 

System. err.println("Erroneous  paragraph  is: "  +  paragraph. getName()); 
System. err.println("Erroneous  line  is: " 

+  paragraph. tagValue("enumerate")): 

System. exit(2); 

} 

if  (enumName. equals(inputStr)) 
break; 

} 

return  enumValue; 

} 

public  void  writeDataFiie(String  fileName) 

{ 

try 

{ 

PrintWriter  pw  =  new  PrintWriter(new  FileWriter(fileName)); 

Iterator  Iterator  =  dataTable.keySet().iterator(): 

while  (iterator.hasNexlO) 

{ 

String  name  =  (String)  iterator.next(); 

Factor  factor  =  (Factor)  dataTable.gel(name); 

int  Q  indexArray  =  new  int  [lnteger.parselnt( 
iniFile.getProfile("config",  "numoffactors"))]; 
for  (int  i  =  0;  i  <  sizeArray[0]:  i++) 

{ 

indexArray[0]  =  i; 

for  (int  j  =  0;  j  <  sizeArray[1];  j++) 

{ 

indexArray[1]  =  j; 
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for  (int  k  =  0;  k  <  sizeArray[2];  k++) 

{ 

indexArray[2]  =  k; 

pw.println(name  +  +  (indexArray[0]  +  1)  + 

+  (indexArray[1]  +  1)  +  +  (indexArray[2]  +  1)  + 

+  factor.get(indexArray)); 

} 

} 

} 

} 

pw.closeO: 

} 

catch  (lOException  ioe) 

{ 

System. err.println(ioe): 

} 

} 

public  Iterator  getKeySetlterator() 

{ 

return  dataTab[e.keySet().iterator(): 

} 

public  Factor  get(String  name) 

{ 

return  (Factor)  dataTable.gel(name): 

} 

public  int  Q  getSizeArrayO 

{ 

return  sizeArray; 

} 

} 
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Factor.iava 


r 

*  Created  on  Aug  28,  2003 
*/ 

package  anova.util; 

r 

*  @author  Purvesh  Khatri 
*1 

import  java.util.ArrayList; 
import  java. util. List; 

public  class  Factor 

{ 

/* 

*  hasChildren  is  only  set  in  the  constructor  and 

*  is  never  modified.  It  is  strictly  a  flag  to  maintain  the  structure 

*  of  the  class  Factor. 

*/ 

private  boolean  hasChildren; 

private  List  childrenList; 

public  Factor(int  Q  sizeArray) 

{ 

childrenList  =  new  ArrayList(sizeArray[0]): 

if  (sizeArray.length  >  1 ) 

{ 

hasChildren  =  true; 

/* 

*  create  new  size  array  for  the  children  of  this  factor  object. 

*  The  new  size  array  removes  the  size  element  from  the  beginning 

*  of  the  size  array. 

V 

int  D  newSizeArray  =  getNewArray(sizeArray): 

//  create  and  add  children  of  this  Factor  instance, 
for  (int  i  =  0;  i  <  sizeArray[0];  i++) 
childrenList.add{i,  new  Factor(newSizeArray)); 

} 

} 

public  Object  get(int  Q  indexArray) 
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{ 

//  If  it  does  not  have  children,  return  the  object  at  given  index, 
if  ( lhasChildren) 

return  childrenList.get(indexArray[0]); 

//Well,  we  have  some  children! 

//  Get  the  child  at  the  given  index  i.e.  indexArray[0]. 

Factor  factor  =  (Factor)  childrenList.get{irdexArray[0]); 

//  Now  reduce  the  indexArray  size  by  one,  i.e., 

II  copy  indices  from  1  to  (length  - 1) 

II  to  a  new  array. 

int  []  newIndexArray  =  getNewArray(indexArray); 

//  Ask  the  child  to  retrieve  the  given  object, 
return  factor.get(newlndexArray); 

} 

private  int  Q  getNewArray(int  Q  array) 

{ 

int  Q  newSizeArray  =  new  int  [array.length  - 1]; 
for  (int  i  =  0;  i  <  newSizeArray.length;  i++j 
newSizeArray[i]  =  array[i  +  1]; 

return  newSizeArray; 

} 

public  void  add(int  □  indexArray,  Double  value) 

{ 

if  ( lhasChildren) 

{ 

childrenList.add(value): 

return; 

} 

Factor  factor  =  (Factor)  childrenList.get(indexArray[0]); 
int  D  newIndexArray  =  getNewArray(indexArray): 
factor.add(newlndexArray,  value); 

} 

public  void  getList(int  Q  array,  List  list) 

{ 

if  (array[0]  ==  -1) 

{ 
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if  ( lhasChildren) 

{ 

for  (int  i  =  0;  i  <  childrenList.size();  i++) 
list.add(childrenList.get(i)); 
return; 

} 

array  =  getNewArray(array): 

for  (int  i  =  0;  i  <  childrenList.size():  i++) 

{ 

Factor  factor  =  (Factor)  childrenList.get(i): 
factor.getList(array,  list); 

} 

return: 

} 

Factor  factor  =  (Factor)  childrenList.get(array[0]); 
int  0  newArray  =  getNewArray(array): 

factor.getList(newArray,  list); 

} 

} 
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IniFile.iava 


package  anova.util; 

import  java. io. File; 
import  java.io.lOException; 
import  java. util. Enumeration; 
import  java. util. Hashtable; 
import  java. util. Set; 

r 

*  this  class  handles  Windows  style  ini-files  of  the  form 

*  [paragraphname] 

*  key1=value1 

*  key2=value2 

* 

*  [another_paragraphname] 

*  keyx  =  vals 

*  tag=foo 

*  etc.  it  will,  however,  work  on  any  platform 
*1 

public  class  IniFile 

{ 

private  String  filename,  path,  fullpath; 
private  Hashtable  paragraphs; 
private  iniParagraph  para; 
private  InputFile  f; 

r 

*  Opens  an  existing  iniFile.  If  it  does  not  exist,  an  empty  one  is  created  and 

*  opened 

* 

*  @param  - 

*  pathname 

*  @param  filename, 

*  usuallyxxx.ini 
*/ 

public  lniFile(String  pathname.  String  fname)  throws  lOException 

{ 

filename  =  fname; 
path  =  pathname; 
if  (path.lengthO  ==  0) 
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fullpath  =  filename; 
else 

fullpath  =  path  +  File. separator  +  filename; 
checkFile(fullpath); 
f  =  new  InputFile(fullpath); 
getParagraphsO; 

} 

private  void  checkFile(String  fname)  throws  lOException 

{ 

File  fl  =  new  Flle{fname): 
if  (  !fl.exists()) 

{ 

OutputFile  outf  =  new  OutputFile(fname); 
outf.close(): 

} 

} 

r 

*  Opens  an  ini  file  in  the  current  directory 

if 

*  @param  - 

*  filename,  usually  xxx.ini 
*1 

public  lniFile{String  fname)  throws  lOException 

{ 

filename  =  fname; 
path  = 

fullpath  =  filename; 
checkFile(fullpath); 
f  =  new  InputFile(filename); 
getParagraphsO; 

} 

private  synchronized  void  getParagraphsO 

{ 

paragraphs  =  new  Hashtable(IO):  //  make  room  for  object  table 
if  ( If.checkErrO) 

{ 

para  =  new  iniParagraph(f); 

if  ( Ipara.errorO) 

{ 

paragraphs. pul(para.getNameO,  para); 

while  ((para.nextLineO  !=  null)  &&  (para.nextLineO-lengthO  >  0)) 

{ 
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para  =  new  iniParagraph(f,  para.nextLine()); 
paragraphs. put(para.getName(),  para); 

} 

} 

f.closeOi 

} 

} 

r* 

*  Returns  an  iniParagraph  instance  for  the  given  paraName. 

* 

*  @param  Paragraph 

*  name. 

*  @return  An  iniParagraph  reference,  if  paraName  is  found,  else  null. 
*! 

public  iniParagraph  getParagraph(String  paraName) 

{ 

return  (iniParagraph)  paragraphs. get{paraName); 

} 

/** 

*  Returns  a  set  of  all  paragraph  names. 

* 

*  ©return  java. util. Set  containing  all  paragraph  sections. 

*/ 

public  Set  getParagraphNamesO 

{ 

return  paragraphs. keySet(); 

} 

private  synchronized  void  putParagraphs(’  throws  lOException 

{ 

OutputFile  f  =  new  OutputFile{fullpath): 

Enumeration  xenum  =  paragraphs. elements(): 
while  (xenum. hasMoreElementsO) 

{ 

iniParagraph  p  =  (iniParagraph)  xenum.nextElement(): 
f.println("["  +  p.getName()  +  "]"); 

Enumeration  e  =  p.getTags(): 
while  (e.hasMoreElements()) 

{ 

String  tag  =  (String)  e.nextElement(): 
f.println(tag.toLowerCase()  +  "="  +  p.tagVa!ue(tag)): 

} 

} 
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f.closeO; 

} 

r 

*  Gets  a  profile  string  from  an  ini  file 

* 

*  @param  para  - 

*  paragraph  name  in  any  case 

*  @param  entry  - 

*  name  of  entry  in  any  case 

*  ©return  String  value  of  this  entry 
*/ 

public  String  getProfile{String  para,  String  entry) 

{ 

iniParagraph  ini  =  null; 

ini  =  (iniParagraph)  paragraphs. get(para.toLowerCase()): 
if  (ini  !=  null)  II  now  look  for  entry  in  that  paragraph 
return  ini.tagValue(entry).trim(); 
else 
return 

} 

r 

*  Gets  a  profile  string  from  an  ini  file  returns  default  value  if  no  such 

*  entry  exists 

* 

*  @param  para  - 

*  paragraph  name  in  any  case  pattern 

*  ©param  entry  - 

*  name  of  entry  in  any  case  pattern 

*  ©param  default  - 

*  default  value  for  parameter 

*  ©return  String  value  of  this  entry 
V 

public  String  getProfile(String  para,  String  entry.  String  defalt) 

{ 

String  value  =  getProfile(para,  entry); 
if  (value. lengthO  <  1) 
value  =  defalt; 
return  value.trim(); 

} 

j** 

*  Puts  a  profile  entry  into  the  ini  file 

* 
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*  @param  para 

*  -paragraph  name 

*  @param  tag  - 

*  name  of  entry 

*  @param  value  - 

*  value  of  entry 
*/ 

public  void  putProfile(String  para,  String  teg,  String  vaiue) 
throws  lOException 

{ 

boolean  found  =  false; 
iniParagraph  p  =  null; 

Enumeration  xenum  =  paragraphs.elements(): 
whiie  (xenum. hasMoreElementsO  &&  Ifound) 

{ 

p  =  (iniParagraph)  xenum. nextElement{); 
found  =  (para.equaisignoreCase(p.getName())); 

} 

if  (found) 

{ 

p.setValue(tag,  value): 

} 

else 

{ 

p  =  new  iniParagraph(para): 
paragraphs. put(p.getName(),  p); 
p.setVaiue(tag.toLowerCase(),  value); 

} 

putParagraphsO:  H  rewrite  entire  ini-file 

} 

} 
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iniParaqraph.lava 

package  anova.util; 

import  java.util.ArrayList; 
import  java. util. Collections; 
import  Java. util. Enumeration: 
import  java. util. Hashtable; 
import  java. util. List; 

/**  A  private  class  used  by  IniFile  7 
public  class  iniParagraph 
{ 

public  List  getSortedTags{) 

{ 

Enumeration  chipEnum  =  getTags(); 

ArrayList  list  =  new  ArrayList(): 

while  (chipEnum. hasMoreElementsO) 

{ 

iist.add((String)  chipEnum. nextElementO); 

} 

Collections.sort(iist): 
return  list; 

} 

private  String  paragraph_name; 

private  InputFile  f; 

private  Hashtable  tags; 

private  String  line; 

private  boolean  error  flag; 

public  iniParagraph{lnputFiie  fl) 

{ 

f  =  fl; 

String  s  =  getNextLine(); 
read_to_nextparagraph(s); 

} 
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public  Enumeration  getTags{) 

{ 

return  tags.keys(): 

) 

public  Hashtabie  getHash{) 

{ 

return  tags; 

} 

public  void  setValue(String  tag,  String  value) 

{ 

tags.put(tag.toLowerCase(),  value);  //  may  replace  old  value 

} 

private  void  read_to_nextparagraph(Strings) 

{ 

error_flag  =  false; 
tags  =  new  Hashtable(); 

//  The  while  loop  looks  for  [paragraph]  line, 
while  ({s  !=  null)  &&  ( !s.startsWith("f’))) 
s  =  getNextLineO; 

//  remove  brackets 
if  (s  !=  null) 

{ 

s  =  s.substring(1 );  //  all  but  first  char 
int  i  =  s.indexOf("]"): 
if  (i  >  0) 

{ 

s  =  s.substring(0,  i); 
store_paragraph(s); 

} 

else 

{ 

error_fiag  =  true; 

} 

} 

else 

{ 

errorflag  =  false; 

} 

} 

public  String  getName() 


56 


{ 

return  paragraph_name.toLowerCase{); 

} 

public  String  tagValue(String  tagname) 

{ 

String  ans; 

ans  =  (String)  tags.get{tagname.toLowerCase()); 
if  (ans  ==  null) 
return 
else 

return  ans; 

} 

private  void  store_paragraph(String  s) 

{ 

String  next; 
inlElement  ini; 
paragraph_name  =  s; 
next  =  getNextLineO; 

while  ((next  !=  null)  &&  ( !next.startsWith("["))) 

{ 

if  (next.startsWithC'#"))//  skip  the  lines  beginning  with  #. 

{ 

next  =  getNextLineO; 
continue; 

} 

ini  =  new  iniElement(next): 
if  ( lini.errorO) 

{ 

II  put  the  property  read  in  the  hash  table. 
tags.put(ini.tagName(),  Ini.valueNameO); 

} 

next  =  getNextLineO: 

} 

if  (tags.sizeO  >  0) 
error_flag  =  false; 

} 

public  iniParagraph(inputFil6  fl,  String  pname) 

{ 

f  =  fl; 

read_to_nextparagraph(pname): 

} 
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public  iniParagraph{String  pname) 

{ 

//  call  this  constructor  to  create  a  new  paragraph 
tags  =  new  Hashtable(): 
paragraph_name  =  pname; 

} 

private  String  getNextLine() 

{ 

line  =  f.readLineO: 

while  ((line  !=  null)  &&  (line,length()  <  1)] 
line  =  f.readLlneQ: 
if  (line  ==  null) 
error_flag  =  true; 
else 

line  =  line.trimO: 
return  line; 

} 

public  boolean  error() 

{ 

return  error_flag: 

} 

public  String  nextLine() 

{ 

return  line; 

} 

} 
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iniElement.iava 

package  anova.util; 

/**  A  private  class  used  by  IniFile  7 
public  class  inlElement 
{ 

private  String  tag,  val; 

private  boolean  error_flag: 

//  parses  one  element  of  an  ini-file 
//  into  a  tag  and  a  value 
public  iniElement(String  s) 

{ 

error_flag  =  false; 
int  i  =  s.indexOf(’'="): 
if  (i  >  0) 

{ 

tag  =  s.substring(0,  i).trim(),toLowerCase(): 
val  =  s.substring(i  +  1); 

} 

else 

errorflag  =  true; 

} 

public  boolean  error() 

{ 

return  error_flag; 

} 

public  String  tagNameQ 

{ 

return  tag; 

} 

public  String  valueNameQ 

{ 

return  val; 

} 

} 
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InputFile.iava 

package  anova.util; 

import  java. io.BufferedReader; 
import  java. io, File; 
import  java.io.FileReader; 
import  java.io.lOException; 

public  class  InputFile 

{ 

private  BufferedReader  f  =  null; 
private  boolean  errflag; 
private  String  s  =  null; 
private  String  filename; 

public  lnputFile(String  fname)  throws  lOException 

{ 

errflag  =  false; 

//  open  file 

fname  =  fname. trlm(): 

f  =  new  BufferedReader(new  FileReader(fname)); 
filename  =  fname; 

} 

public  boolean  checkErr() 

{ 

return  errflag; 

} 

public  String  read() 

{ 

//  read  a  single  field  up  to  a  comma  or  end  of  line 
String  ret  = 
if  (s  ==  null) 

{ //  if  no  data  in  string 
s  =  readLine();  //  read  next  line 

} 

if  (s  !=  null) 

{ //  if  there  is  data 
s.trim();  //  trim  off  blanks 
int  i  =  s.indexOf(",");  //  find  next  comma 
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if(i<=0) 

{ 

ret  =  s.trimO;  //  if  no  commas  go  to  end  of  line 
s  =  null;  II  and  null  out  stored  string 

} 

else 

{ 

ret  =  s,substring(0,  i).trim();  //  return  left  of  comma 
s  =  s.substring(i  +  1 ):  //  save  right  of  comma 

} 

} 

else 

ret  =  null; 

return  ret;  II  return  string 

} 

public  String  readLine() 

{ 

//  read  in  a  line  from  the  file 
s  =  null; 
try 
{ 

s  =  f.readLineO:  II  could  throw  error 

} 

catch  (lOException  e) 

{ 

errflag  =  true; 

System.out.printlnC'FlIe  read  error"); 

} 

return  s; 

} 

public  void  close{) 

{ 

try 

{ 

f.closeO;  II  close  file 

} 

catch  (lOException  e) 

{ 

System, outprintlnC'File  close  error"); 
errflag  =  true; 

} 

} 

public  String  root() 
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{ 

int  i  =  filename. lastlndexOf(File.separatorChar): 
String  s  =  filename.substring{i  +  1 ); 
i  =  s.indexOf("."): 
return  s.substring{0,  i); 

} 

} 
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OutputFiie.iava 

package  anova.util; 

import  java.io.BufferedWriter; 
import  java.io.FileWriter; 
import  java.io.lOException; 
import  java.io.PrintWriter; 

//  Output  file  class 
pubiic  ciass  OutputFile 
{ 

private  BufferedWriter  f; 
private  PrintWriter  p; 
private  int  tabcolumn; 
private  int  width; 

r 

*  Create  output  file 

* 

*  @param  complete 

*  path  and  name  of  file 
*/ 

public  OutputFile(String  filename)  throws  lOException 

{ 

tabcolumn  =  0; 
width  =  0; 

fiiename  =  filename.trim(): 

System. out.println(filename); 
f  =  new  BufferedWriter(new  FileWriter(filename)): 
p  =  new  PrintWriter(f): 

} 

/**  insert  a  tab  into  the  file  */ 
public  void  tab() 

{ 

p.print("\t"): 

} 

jit* 

*  insert  spaces  into  the  file 

* 

*  @param  the 
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number  of  spaces  to  insert 


it 

*/ 

public  String  space(int  n) 

{ 

StringBuffer  sb  =  new  StringBuffer(n): 

//  put  spaces  into  string  buffer 
for  (int  i  =  0;  i  <  n;  i++) 

{ 

sb.insert(i, ' '); 

} 

return  sb.toStringO; 

} 

r 

*  tab  over  to  column 

* 

*  @param  the 

*  column  to  tab  over  to  Inserts  spaces  to  move  up  to  tab  column  or 

*  starts  a  new  line 
*/ 

public  void  tab(int  tb) 

{ 

if  (tb  >  tabcolumn) 

{ 

print(space(tb  -  tabcolumn)); 

} 

else 

printInC'"); 

} 

public  void  println(String  s) 

{ 

p.println(s); 
tabcolumn  =  0; 

} 

public  void  print!n(int  i) 

{ 

p.println(i): 
tabcolumn  =  0; 

} 

public  void  println(double  d) 

{ 

tabcolumn  =  0; 
p.println(d); 
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} 


public  void  print(String  s) 

{ 

p.print{s): 

tabcolumn  +=  s.length(): 

} 

public  void  print(int  i) 

{ 

String  s  =  new  lnteger(i).toString(): 
if  (s.length()<  width) 
print(space(width  -  s.length())): 
print(s); 

} 

public  void  print(float  f) 

{ 

String  s  =  new  Fioat{f).toString(); 
print(s): 

} 

public  void  print(double  d) 

{ 

String  s  =  new  Double(d).toString(): 
print(s): 

} 

public  void  close() 

{ 

p.closeO; 

} 

public  void  finalizeQ 

{ 

close(); 

} 

} 
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APPENDIX  C 


The  following  files  can  be  used  to  test  an  installation  of  the  ANOVA  Java  program  to 
insure  that  it  is  operating  properly. 

TEST  CONFIGURATION  FILE 

[config] 

#  Number  of  factors  to  consider  in  Anova. 

NumOfFactors  =  3 

#  The  following  must  be  enclosed  in  double  quotes. 

ColumnSeparator  = 

DataFile  =  sampleinput.txt 

TwoWayOutputFile  =  SampleTwoWayAnovaOutput.txt 
OneWayOutputFile  =  SampleOneWayAnovaOutput.txt 
HeadersForOneWayAnova  =  Factor2 

#  Now  describe  each  factor  individually. 

#  A  factor  paragraph  has  the  format  Factor#,  where  #  is  substituted  by  an  integer. 
[Factorl] 

Name  =  sample 
Columns  =  2 
MaxNumOfIndex  =  6 

[Factor2] 

Name  =  condition 
Columns  =  3 
MaxNumOfIndex  =  3 

Enumerate  =  Heat  1  Cold  2  Control  3 

[Factor3] 

Name  =  time 
Columns  =  4 
MaxNumOfIndex  =  2 
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TEST  INPUT  FILE 


AFFX-BioB-5_at,1  .Control,  1 ,470.90 
AFFX-BioB-5_at,1,Control,2,455.50 
AFFX-BioB-5_at,1  .Heat,  1 ,486.70 
AFFX-BioB-5_at,1,Heat,2,464.40 
AFFX-BioB-5_at,1  .Cold, 1 ,378.40 
AFFX-BioB-5_at,1  ,Cold.2, 380.00 
AFFX-BioB-5_at,2, Control, 1,370.40 
AFFX-BioB-5_at, 2, Control,2, 399.30 
AFFX-BioB-5_at,2,Heat,1 ,41 8.20 
AFFX-BioB-5_at,2,Heat,2,454.70 
AFFX-BioB-5_at,2,Cold,1 ,395.30 
AFFX-BioB-5_at,2,Cold,2, 425.40 
AFFX-BioB-5_at,3,Control, 1,393.40 
AFFX-BioB-5_at,3,Control,2, 392.70 
AFFX-BioB-5_at,3,Heat,1 ,386.60 
AFFX-BioB-5_at,3,Heat,2,428.60 
AFFX-BioB-5_at,3,Cold,1 ,282.90 
AFFX-BioB-5_at,3,Cold,2,331 .80 
AFFX-BioB-5_at, 4, Control, 1,347.80 
AFFX-BioB-5_at, 4, Control, 2,389.80 
AFFX-BioB-5_at, 4, Heat, 1,393.80 
AFFX-BloB-5_at, 4, Heat,2, 443.60 
AFFX-BioB-5_at,4,Cold,1 ,354,20 
AFFX-BioB-5_at,4,Cold,2,432.30 
AFFX-BioB-5_at.5, Control, 1,339.50 
AFFX-BioB-5_at,5,Control,2,373.10 
AFFX-BioB-5_at,5,Heat,1 ,402.40 
AFFX-BioB-5_at,5,Heat,2, 524.00 
AFFX-BioB-5_at,5,Cold,1 ,346.1 0 
AFFX-BioB-5_at.5.Cold, 2,360.60 
AFFX-BioB-5_at, 6, Control,  1,354.50 
AFFX-BioB-5_at, 6, Control,2, 314.40 
AFFX-BioB-5_at,6,Heat,1 ,376.80 
AFFX-BioB-5_at, 6, Heat, 2, 1097.10 
AFFX-BioB-5_at, 6, Cold, 1,484.50 
AFFX-BioB-5_at,6,Cold,2, 419.20 
AFFX-BioB-M_at,1  .Control, 1 ,776.20 
AFFX-BioB-M_at,1, Control, 2, 963.80 
AFFX-BioB-M_at.1  ,Heat,  1 ,734.30 
AFFX-BioB-M_at,1.Heat,2,700.20 
AFFX-BioB-M_at,  1  ,Cold.1 ,630.00 
AFFX-BioB-M_at,  1  ,Coid.2,641 .00 
AFFX-BioB-M_at, 2, Control, 1,649.10 
AFFX-BioB-M_at, 2, Control, 2, 787.90 
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AFFX-BioB-M_at, 2, Heat, 1 ,627.60 
AFFX-BioB-M_at,2,Heat,2,528.20 
AFFX-BioB-M_at,2,Cold,1 ,751 .40 
AFFX-BioB-M_at.2,Cold,2,501.10 
AFFX-BioB-M_at,3,Control,1 ,589.00 
AFFX-BioB-M_at,3,Control, 2,655.50 
AFFX-BioB-M_at,3,Heat,1 ,720.30 
AFFX-BioB-M_at,3,Heat,2,534.00 
AFFX-BioB-M_at,3,Cold,1 ,674.80 
AFFX-BioB-M_at,3,Cold,2, 552,00 
AFFX-BioB-M_at,4, Control, 1,400. 70 
AFFX-BioB-M_at,4,Control,2,479.40 
AFFX-BioB-M_at,4,Heat,1 ,520.60 
AFFX-BioB-M_at,4,Heat,2,61 6.80 
AFFX-BioB-M_at,4,Cold,  1 ,423.60 
AFFX-BioB-M_at,4,Cold,2, 619.40 
AFFX-BioB-M_at,5, Control, 1 ,431 .50 
AFFX-BioB-M_at,5,Control,2, 530.00 
AFFX-BioB-M_at,5,Heat,1 ,578.40 
AFFX-BioB-M_at,5,Heat,2,661.10 
AFFX-BioB-M_at,5,Cold,1 ,450.80 
AFFX-BioB-M_at,5, Cold, 2,462.60 
AFFX-BioB-M_at,6,Control,  1 ,535.40 
AFFX-BioB-M_at.6.Control.2.389.90 
AFFX-BioB-M_at,6,Heat,1 ,529.90 
AFFX-BioB-M_at,6,Heat,2,1789.60 
AFFX-BioB-M_at,6,Cold, 1,558.20 
AFFX-BioB-M_at,6,Cold,2,656.60 
AFFX-BioB-3_at,1  ,Control,1 ,405.00 
AFFX-BioB-3_at,1  ,Control,2,422.1 0 
AFFX-BioB-3_at,1  ,Heat,1 ,475.60 
AFFX-BioB-3_at,1  ,Heat,2,537.20 
AFFX-BioB-3_at,1  ,Cold,1 ,396.50 
AFFX-BioB-3_at,1  ,Cold,2,444.00 
AFFX-BioB-3_at,2,Control, 1,202. 10 
AFFX-BioB-3_at,2,Control, 2,436,90 
AFFX-BioB-3_at,2,Heat,1, 364.70 
AFFX-BioB-3_at,2.Heat,2,365.80 
AFFX-BioB-3_at,2,Cold,1 ,393.1 0 
AFFX-BioB-3_at,2,Cold,2,274.30 
AFFX-BioB-3_at,3,  Control,  1 ,308.80 
AFFX-BloB-3_at,3,Control, 2,425.60 
AFFX-BloB-3_at, 3, Heat,1, 516.90 
AFFX-BioB-3_at,3,Heat,2, 368.60 
AFFX-BioB-3_at,3,Cold,1 ,384.70 
AFFX-BloB-3_at,3,Cold,2,320.50 
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AFFX-BioB-3_at, 4, Control, 1,249.00 
AFFX-BioB-3_at, 4, Control,2, 265.60 
AFFX-BioB-3_at,4,Heat,1 ,272.60 
AFFX-BioB-3_at,4,Heat,2,335.40 
AFFX-BioB-3_at, 4, Cold, 1 ,255.00 
AFFX-BioB-3_at,4,Cold, 2, 290.40 
AFFX-BioB-3_at,5,Control,  1 ,293.80 
AFFX-BioB-3_at,5, Control, 2,252.30 
AFFX-BioB-3_at, 5, Heat,1, 304.50 
AFFX-BioB-3_at,5, Heat, 2, 400.50 
AFFX-BioB-3_at,5,Cold,1 ,260.80 
AFFX-BioB-3_at,5,Cold, 2, 266.20 
AFFX-BioB-3_at,6,Control, 1,286.30 
AFFX-BioB-3_at,6,Control,2,278.70 
AFFX-BioB-3_at,6,Heat, 1,303.80 
AFFX-BioB-3_at,6,Heat,2, 1142.50 
AFFX-BioB-3_at,6,Cold,1 ,307.60 
AFFX-BioB-3_at,6,Cold,2,31 2.00 
AFFX-BioC-5_at,1, Control, 1,1 036.40 
AFFX-BioC-5_at,1  .Control, 2,1 1 05.70 
AFFX-BloC-5_at,1  ,Heat,1 ,1 051 .80 
AFFX-BloC-5_at,1,Heat,2,982.60 
AFFX-BioC-5_at,1, Cold, 1,1066.50 
AFFX-BioC-5_at,1 , Cold, 2,1 1 70.70 
AFFX-BioC-5_at, 2, ControU, 907.40 
AFFX-BioC-5_at,2, Control, 2, 986.50 
AFFX-BioC-5_at,2,Heat,1 ,1 144.60 
AFFX-BioC-5_at,2,Heat,2,1067.00 
AFFX-BioC-5_at,2,Cold,1 ,1 1 1 3.90 
AFFX-BioC-5_at,2,Cold, 2, 909.20 
AFFX-BioC-5_at,3,Control,  1 ,91 8.30 
AFFX-BioC-5_at,3,Control,2,954.30 
AFFX-BioC-5_at,3, Heat, 1,1032.80 
AFFX-BloC-5_at,3,Heat,2,864.50 
AFFX-BioC-5_at,3,Cold,  1 ,950.50 
AFFX-BioC-5_at,3,Cold,2.803.30 
AFFX-BioC-5_at,4,Control, 1,872.50 
AFFX-BioC-5_at,4,Control, 2, 923.40 
AFFX-BioC-5_at,4,Heat,1 ,939.60 
AFFX-BioC-5_at,4,Heat,2,1 033.90 
AFFX-BioC-5_at,4,Cold,1 ,857.80 
AFFX-BioC-5_at,4,Cold,2,986.20 
AFFX-BioC-5_at, 5, ControU, 1067.30 
AFFX-BioC-5_at,5,Control,2,1071.00 
AFFX-BioC-5_at,5,Heat,1 ,1 122.30 
AFFX-BioC-5_at,5,Heat,2, 1226.00 
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AFFX-BioC-5_at,5,Cold,1 ,91 1 .80 
AFFX-BioC-5_at,5,Cold.2,1 01 0.70 
AFFX-BioC-5_at, 6, Control, 1,936.10 
AFFX-BioC-5_at, 6, Control,2, 974.40 
AFFX-BioC-5_at,6,Heat, 1,943.00 
AFFX-BioC-5_at, 6, Heat, 2,2778.90 
AFFX-BioC-5_at,6,Cold,1 ,1 150.00 
AFFX-BioC-5_at, 6, Cold, 2,1 222.30 
AFFX-BioC-3_at,1 , Control,  1 ,989. 1 0 
AFFX-BioC-3_at,1, Control, 2, 870.20 
AFFX-BioC-3_at,1,Heat,1. 910.80 
AFFX-BioC-3_at,1 , Heat, 2, 879.90 
AFFX-BioC-3_at,1, Cold, 1,820.10 
AFFX-BioC-3_at,1 , Cold, 2,81 6.60 
AFFX-BloC-3_at,2,Control, 1,668.00 
AFFX-BioC-3_at,2,Control,2, 749.70 
AFFX-BioC-3_at,2,Heat,1 ,864.30 
AFFX-BioC-3_at,2,Heat,2, 759.30 
AFFX-BioC-3_at, 2, Cold,1, 824.00 
AFFX-BioC-3_at,2,Cold,2,731 .30 
AFFX-BioC-3_at,3,Control, 1,575.10 
AFFX-BioC-3_at, 3, Control, 2,741 .90 
AFFX-BioC-3_at,3,Heat,1 ,801 .40 
AFFX-BioC-3_at,3,Heat,2.699.80 
AFFX-BioC-3_at,3,Cold,1 ,758.40 
AFFX-BioC-3_at,3,Cold,2,641 .60 
AFFX-BioC-3_at, 4, Control,  1 ,731 .50 
AFFX-BioC-3_at,4,Control,2, 734.40 
AFFX-BioC-3_at,4,Heat,  1 .739.00 
AFFX-BioC-3_at,4,Heat, 2,788.30 
AFFX-BtoC-3_at,4,Cold, 1,696.30 
AFFX-BioC-3_at,4,Cold,2,771 .10 
AFFX-BloC-3_at,5,Control,1,851.50 
AFFX-BioC-3_at,5,Control,2,838.10 
AFFX-BioC-3_at,5,Heat,1 ,804.50 
AFFX-BloC-3_at,5.Heat, 2, 1014.00 
AFFX-BioC-3_at, 5, Cold,1, 692.00 
AFFX-BioC-3_at, 5, Cold, 2,810.70 
AFFX-BioC-3_at, 6, Control, 1,752.20 
AFFX-BioC-3_at, 6, Control, 2, 827.40 
AFFX-BioC-3_at,6,Heat,1 ,753.40 
AFFX-BioC-3_at,6,Heat,2,281 0.00 
AFFX-BioC-3_at,6,Cold,  1 ,843.90 
AFFX-BioC-3_at, 6, Cold, 2,770.60 
AFFX-BioDn-5_at,1 , Control, 1 ,101 3.50 
AFFX-BioDn-5_at,1  ,Control,2,1 263.80 
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AFFX-BioDn-5_at.1  ,Heat,1 ,1 331 .20 
AFFX-BioDn-5_at,1, Heat, 2, 1157.70 
AFFX-BioDn-5_at,1, Cold, 1,1008.00 
AFFX-BioDn-5_at,1  .Cold, 2,935.00 
AFFX-BioDn-5_at,2,Control,1, 867.00 
AFFX-BioDn-5_at,2,Control,2,1007.60 
AFFX-BioDn-5_at,2,Heat,1, 1132.10 
AFFX-BioDn-5_at,2,Heat, 2,1013.60 
AFFX-BioDn-5_at,2,Cold,1 ,927.50 
AFFX-BioDn-5_at,2,Cold,2,966.40 
AFFX-BioDn-5_at,3,Control.  1 .885.50 
AFFX-BioDn-5_at,3,Control,2,881.40 
AFFX-BioDn-5_at,3,Heat, 1,812.20 
AFFX-BioDn-5_at,3,Heat, 2, 950.00 
AFFX-BioDn-5_at, 3, Cold,1 ,885.90 
AFFX-BioDn-5_at,3,Cold, 2.728.20 
AFFX-BioDn-5_at,4.Control, 1,949.10 
AFFX-BioDn-5_at,4,Control,2,1097.60 
AFFX-BioDn-5_at,4,Heat,1 ,1 040.70 
AFFX-BioDn-5_at,4,Heat,2,1344.60 
AFFX-BioDn-5_at,4,Cold,1 ,929.30 
AFFX-BioDn-5_at,4,Cold,2,1 1 35.10 
AFFX-BioDn-5_at,5, Control,  1 , 1 247.80 
AFFX-BioDn-5_at,5,Control,2, 1 068.40 
AFFX-BioDn-5_at,5,Heat, 1,1262.20 
AFFX-BioDn-5_at,5,Heat,2, 1570.90 
AFFX-BioDn-5_at,5,Cold,1 .1 172.80 
AFFX-BioDn-5_at,5,Cold, 2, 1060.80 
AFFX-BioDn-5_at,6.Control,1 ,1 053.60 
AFFX-BioDn-5_at,6,Control,2,1086.10 
AFFX-BioDn-5_at,6,Heat,1 ,1 146.90 
AFFX-BioDn-5_at,6,Heat,2,3057.50 
AFFX-BioDn-5_at,6,Cold,1 ,1 1 86.40 
AFFX-BioDn-5_at,6,Cold,2,1 1 86.20 
AFFX-BioDn-3_at,  1  .Control,  1 ,7207.30 
AFFX-BioDn-3_at,1,Control,2, 9109.50 
AFFX-BioDn-3_at,1  .Heat,1 ,7460.60 
AFFX-BioDn-3_at,1,Heat,2,8004.00 
AFFX-BioDn-3_at,1  ,Cold,1 ,6797.00 
AFFX-BioDn-3_at,1, Cold, 2,7549.40 
AFFX-BioDn-3_at.2,Control,1, 6300.50 
AFFX-BioDn-3_at,2, Control, 2, 6465.00 
AFFX-BioDn-3_at, 2, Heat, 1.7307.40 
AFFX-BioDn-3_at,2, Heat, 2, 6776.60 
AFFX-BioDn-3_at,2,Cold,1, 6472.10 
AFFX-BioDn-3  at,2,Cold,2, 6229.00 
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AFFX-BioDn-3_at, 3, Control,  1 ,6437.40 
AFFX-BioDn-3_at,3,Control,2,6234.80 
AFFX-BioDn-3_at,3,Heat, 1,6924.40 
AFFX-BioDn-3_at,3,Heat,2,6761.00 
AFFX-BioDn-3_at, 3, Cold, 1,6080.10 
AFFX-BioDn-3_at,3,Cold,2,6227.30 
AFFX-BioDn-3_at,4,Control.1,6211.50 
AFFX-BioDn-3_at, 4, Control, 2, 5820.80 
AFFX-BioDn-3_at,4.Heat, 1,6362.70 
AFFX-BioDn-3_at,4,Heat,2,8059.10 
AFFX-BioDn-3_at, 4, Cold, 1,6298.40 
AFFX-BioDn-3_at,4,Cold, 2, 6755.40 
AFFX-BioDn-3_at,5, Control, 1,6650.00 
AFFX-BioDn-3_at, 5, Control, 2, 6296.00 
AFFX-BioDn-3_at, 5, Heat,1, 7538.70 
AFFX-BioDn-3_at, 5, Heat,2, 9636.70 
AFFX-BioDn-3_at, 5, Cold, 1,6252.10 
AFFX-BioDn-3_at.5.Cold,2,6886.40 
AFFX-BioDn-3_at,6,Control,1, 7114.50 
AFFX-BioDn-3_at,6,Control,2,6865.70 
AFFX-BioDn-3_at,6,Heat,1 ,81 75.40 
AFFX-BioDn-3_at,6,Heat, 2,22718.70 
AFFX-BioDn-3_at,6,Cold,1,8222.10 
AFFX-BioDn-3_at,6,Cold,2.7978.20 
AFFX-CreX-5_at,1 , Control,!  ,10349.20 
AFFX-CreX-5_at,1,Control,2, 10896.80 
AFFX-CreX-5_at,1,Heat,1, 10961. 60 
AFFX-CreX-5_at,1,Heat,2,10869.30 
AFFX-Cr6X-5_at,1, Cold, 1,10200.50 
AFFX-CreX-5_at,1, Cold, 2, 9888.20 
AFFX-CreX-5_at,2, Control,!  ,821 8.50 
AFFX-CreX-5_at,2,Control,2,1 1 730.00 
AFFX-CreX-5_at,2,Heat,1 ,1 0366.00 
AFFX-CreX-5_at,2,Heat, 2, 10479.70 
AFFX-CreX-5_at,2,Cold, 1,10053.30 
AFFX-CreX-5_at,2,Cold,2,10261.80 
AFFX-CreX-5_at,3,Control,  1 ,9099.40 
AFFX-CreX-5_at,3,Contro!, 2, 10479.70 
AFFX-CreX-5_at,3,Heat,  1 , 1 0265.50 
AFFX-CreX-5_at,3,Heat,2, 9779.50 
AFFX-CreX-5_at, 3, Cold, 1,8689.50 
AFFX-CreX-5_at,3,Cold.2, 10807.20 
AFFX-CreX-5_at, 4, Control,!, 8572.30 
AFFX-CreX-5_at,4,Control,2,8808.80 
AFFX-CreX-5_at,4,Heat,1 ,9909.80 
AFFX-CreX-5_at,4,Heat,2,1 11 75.80 
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AFFX-CreX-5_at, 4, Cold,  1 ,8380.30 
AFFX-CreX-5_at,4,Cold,2, 1 0343.70 
AFFX-CreX-5_at, 5, Control, 1,10793.60 
AFFX-CreX-5_at,5,Control.2, 9764.40 
AFFX-CreX-5_at, 5, Heat,1, 11415.80 
AFFX-CreX-5_at,5,Heat,2, 14937.50 
AFFX-CreX-5_at,5,Cold,1 ,9218.20 
AFFX-CreX-5_at,5, Cold, 2, 11477.30 
AFFX-CreX-5_at,6,Control,1, 10291. 30 
AFFX-CreX-5_at,6, Control, 2,9882.50 
AFFX-CreX-5_at,6,Heat,1 ,1 0594.00 
AFFX-CreX-5_at, 6, Heat, 2,31006.60 
AFFX-CreX-5_at,6,Cold,1, 10829.10 
AFFX-CreX-5_at, 6, Cold, 2, 11607.30 
AFFX-CreX-3_at,1, Control, 1,14407.10 
AFFX-CreX-3_at,1 , Control, 2,1 6371 .70 
AFFX-CreX-3_at,1 , Heat, 1 ,14402.70 
AFFX-CreX-3_at,1,Heat,2, 15447.60 
AFFX-CreX-3_at,1, Cold, 1,14135.00 
AFFX-CreX-3_at,1,Cold,2,15241.30 
AFFX-CreX-3_at, 2, Control,  1,11 974.20 
AFFX-CreX-3_at,2,Control,2,15715.90 
AFFX-CreX-3_at, 2, Heat, 1,15349.30 
AFFX-CreX-3_at, 2, Heat,2, 16210.30 
AFFX-CreX-3_at, 2, Cold, 1,14645.00 
AFFX-CreX-3_at,2,Cold,2, 14767.40 
AFFX-CreX-3_at,3,Control, 1,12924.20 
AFFX-CreX-3_at.3, Control, 2,14047.80 
AFFX-CreX-3_at,3,Heat,1 ,1 5006.80 
AFFX-CreX-3_at,3,Heat, 2, 14172.70 
AFFX-CreX-3_at,3,Cold,1, 12170.40 
AFFX-CreX-3_at,3,Cold, 2,14946.40 
AFFX-CreX-3_at,4,Control.1 ,1 1 1 79.20 
AFFX-CreX-3_at, 4, Control, 2,1 1 960.40 
AFFX-CreX-3_at,4,Heat, 1,12972.60 
AFFX-CreX-3_at,4,Heat, 2, 14237.70 
AFFX-CreX-3_at, 4, Cold.1. 11492.30 
AFFX-CreX-3_at,4,Cold,2,13275.60 
AFFX-CreX-3_at, 5, Control, 1,15797.40 
AFFX-CreX-3_at,5,Control,2,1 541 0.30 
AFFX-CreX-3_at,5,Heat,1 ,1 6493.20 
AFFX-CreX-3_at,5, Heat, 2, 20345.00 
AFFX-CreX-3_at,5,Cold,1 ,1 2950.00 
AFFX-CreX-3_at,5,Cold, 2,16059.10 
AFFX-CreX-3_at,6,Control,  1 ,1 3583.80 
AFFX-CreX-3_at,6,Control,2, 1 4567.80 
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AFFX-CreX-3_at, 6, Heat, 1,14343.30 
AFFX-CreX-3_at,6,Heat,2,41 726. 1 0 
AFFX-CreX-3_at.6,Cold,1 ,1 5294.30 
AFFX-CreX-3_at, 6, Cold, 2, 14518.30 
AFFX-DapX-5_at,1 , Control, 1 ,14.30 
AFFX-DapX-5_at,1  ,Control,2,61 .50 
AFFX-DapX-5_at,1 , Heat, 1 ,24.70 
AFFX-DapX-5_at,1 , Heat, 2,72.30 
AFFX-DapX-5_at,1  ,Cold,1 ,45.50 
AFFX-DapX-5_at,1, Cold, 2,90.00 
AFFX-DapX-5_at, 2, Control,  1 ,93.60 
AFFX-DapX-5_at,2,Control,2,74.40 
AFFX-DapX-5_at,2,Heat,1,37.80 
AFFX-DapX-5_at,2,Heat,2,86.20 
AFFX-DapX-5_at,2,Cold,1,56.10 
AFFX-DapX-5_at,2,Cold,2,33. 1 0 
AFFX-DapX-5_at,3,Control,1 ,9.90 
AFFX-DapX-5_at,3,Control, 2,42.00 
AFFX-DapX-5_at,3,Heat,1 ,39.50 
AFFX-DapX-5_at,3,Heat, 2,23.50 
AFFX-DapX-5_at,3,Cold, 1,3.90 
AFFX-DapX-5_at,3,Cold,2,53.80 
AFFX-DapX-5_at,4,Control, 1,80.00 
AFFX-DapX-5_at.4,Control.2,8.00 
AFFX-DapX-5_at,4,Heat,1, 56.50 
AFFX-DapX-5_al,4,Heat,2,41.10 
AFFX-DapX-5_at,4,Cold,1 ,76.30 
AFFX-DapX-5_at,4,Cold, 2,65.10 
AFFX-DapX-5_at, 5, Control, 1,36. 00 
AFFX-DapX-5_at,5,Control,2,26.10 
AFFX-DapX-5_at,5,Heat,1,9.80 
AFFX-DapX-5_at,5,Heat,2,36.50 
AFFX-DapX-5_at,5,Cold,1 ,27.60 
AFFX-DapX-5_at,5,Cold,2,12.80 
AFFX-DapX-5_at,6,Control,1, 30.40 
AFFX-DapX-5_at,6,Control, 2,14.50 
AFFX-DapX-5_at,6,Heat,1 ,38.00 
AFFX-DapX-5_at,6,H6at, 2,1 5.80 
AFFX-DapX-5_at,6,Cold,1. 49.30 
AFFX-DapX-5_at,6,Cold,2,7.80 
AFFX-DapX-M_at,1 , Control, 1 ,66.30 
AFFX-DapX-M_at,  1 ,  Control  ,2 ,116.80 
AFFX-DapX-M_at,1  ,Heat,1 ,77.60 
AFFX-DapX-M_at,1,Heat,2,45.80 
AFFX-DapX-M_at,  1 ,  Cold  ,1,113.50 
AFFX-DapX-M_at,1,Cold,2, 65.40 


AFFX-DapX-M_at, 2, Control,!  ,1 33.40 
AFFX-DapX-M_at, 2, Control, 2, 89.70 
AFFX-DapX-M_at,2,Heat,1 ,56.30 
AFFX-DapX-M_at,2,Heat,2,53.90 
AFFX-DapX-M_at,2,Cold.1, 60.10 
AFFX-DapX-M_at,2,Cold,2, 29.30 
AFFX-DapX-M_at,3,Control, 1,56.70 
AFFX-DapX-M_at,3,Control,2,142.60 
AFFX-DapX-M_at,3,Heat,1 ,54.00 
AFFX-DapX-M_at,3,Heat,2,86.70 
AFFX-DapX-M_at, 3, Cold, 1.1 13.40 
AFFX-DapX-M_at.3,Cold,2, 10.40 
AFFX-DapX-M_at, 4, Control,  1 ,92.30 
AFFX-DapX-M_at.4,Control.2.46.10 
AFFX-DapX-M_at,4,Heat,1 ,93.00 
AFFX-DapX-M_at.4,Heat,2,81 .00 
AFFX-DapX-M_at,4,Cold,1 ,1 09.70 
AFFX-DapX-M_at,4,Cold,2,35.00 
AFFX-DapX-M_at,5,Control,  1 ,66.20 
AFFX-DapX-M_at.5,Control.2.66.30 
AFFX-DapX-M_at,5,Heat,1 .9.80 
AFFX-DapX-M_at,5, Heat, 2,72.30 
AFFX-DapX-M_at,5.Cold, 1,1 02.10 
AFFX-DapX-M_at,5,Cold,2,1 01 .50 
AFFX-DapX-M_at,6,Control, 1,78.30 
AFFX-DapX-M_at,6,Control,2,57.90 
AFFX-DapX-M_at,6,Heat,1, 127.80 
AFFX-DapX-M_at,6,Heat,2,221 .80 
AFFX-DapX-M_at,6,Cold,1 ,94.20 
AFFX-DapX-M_at,6.Cold, 2,99.10 
AFFX-DapX-3_at,1  .Control,!  ,35.50 
AFFX-DapX-3_at,1  ,Control,2,1 1 .00 
AFFX-DapX-3_at,1 , Heat,!  ,7.80 
AFFX-DapX-3_at,  1 ,  Heat,  2, 1 4. 1 0 
AFFX-DapX-3_at,  1  .Cold  ,1,11.20 
AFFX-DapX-3_at,1 , Cold, 2,1 1 .00 
AFFX-DapX-3_at.2.Control,1 ,1 1 .80 
AFFX-DapX-3_at,2,Control,2,17.80 
AFFX-DapX-3_at,2,Heat,1,9.00 
AFFX-DapX-3_at,2,Heat,2,12.30 
AFFX-DapX-3_at,2,Cold, 1,8.50 
AFFX-DapX-3_at,2,Cold,2,20.40 
AFFX-DapX-3_at,3,Control,1 ,1 5.50 
AFFX-DapX-3_at,3.Control,2,12.40 
AFFX-DapX-3_at,3, Heat,!, 9.50 
AFFX-DapX-3_at,3,Heat,2,7.30 
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AFFX-DapX-3_at,3, Cold, 1 ,5.60 
AFFX-DapX-3_at,3,Coid,2,14.90 
AFFX-DapX-3_at, 4, Control,  1 , 1 1 .60 
AFFX-DapX-3_at,4,Contro!,2,7.60 
AFFX-DapX-3_at,4,Heat,  1 ,9.40 
AFFX-DapX-3_at,4,Heat,2,1 1 .40 
AFFX-DapX-3_at,4,Cold, 1,1 2.10 
AFFX-DapX-3_at,4,Cold,2,10.80 
AFFX-DapX-3_at, 5, Control, 1,8. 20 
AFFX-DapX-3_at, 5, Control, 2, 8.70 
AFFX-DapX-3_at, 5, Heat, 1,6.60 
AFFX-DapX-3_at,5,Heat,2,1 1.10 
AFFX-DapX-3_at, 5, Cold, 1,10.90 
AFFX-DapX-3_at,5,Cold,2,12.90 
AFFX-DapX-3_at,6,Control,1,12.70 
AFFX-DapX-3_at,6,Control,2,7.60 
AFFX-DapX-3_at,6,Heat,1 ,5.80 
AFFX-DapX-3_at,6,Heat,2,18.70 
AFFX-DapX-3_at,6,Cold,  1 ,1 3.40 
AFFX-DapX-3_at,6,Cold,2,9.00 
AFFX-LysX-5_at,1 , Control, 1 ,5.30 
AFFX-LysX-5_at,1, Control, 2, 29.70 
AFFX-LysX-5_at,1  ,Heat,1 ,6.40 
AFFX-LysX-5_at,1 , Heat, 2, 6.30 
AFFX-LysX-5_at,1  .Cold,  1 ,28.90 
AFFX-LysX-5_at,1  .Cold, 2,57.20 
AFFX-LysX-5_at,2,Control, 1,44.10 
AFFX-LysX-5_at,2,Control,2,32.10 
AFFX-LysX-5_at,2,Heat,1 ,8.70 
AFFX-LysX-5_at,2,Heat,2,2.50 
AFFX-LysX-5_at,2,Cold,  1 ,44.30 
AFFX-LysX-5_at,2,Cold,2,4.20 
AFFX-LysX-5_at,3,Control, 1,14.20 
AFFX-LysX-5_at,3,Control, 2,30.70 
AFFX-LysX-5_at,3,Heat,1 ,25.70 
AFFX-LysX-5_at,3,Heat,2,12.30 
AFFX-LysX-5_at,3,Cold,  1 ,20.40 
AFFX-LysX-5_at,3,Cold,2,9.50 
AFFX-LysX-5_at,4, Control, 1 ,7.00 
AFFX-LysX-5_at,4,Control,2,8.00 
AFFX-LysX-5_at,4,Heat,1 ,22.60 
AFFX-LysX-5_at,4,Heat,2,6.40 
AFFX-LysX-5_at,4,Cold,1 ,5.70 
AFFX-LysX-5_at,4,Cold,2,16.30 
AFFX-LysX-5_at,5,Control,1, 20.50 
AFFX-LysX-5_at,5,Control,2,21.10 
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AFFX-LysX-5_at,5,  Heat,  1,12.40 
AFFX-LysX-5_at,5.Heat,2,9.80 
AFFX-LysX-5_at.5,Cold,1,6.10 
AFFX-LysX-5_at,5,Cold,2,8.60 
AFFX-LysX-5_at,6, Control, 1 ,1 1 .00 
AFFX-LysX-5_at,6,Control,2,6.80 
AFFX-LysX-5_at,6,Heat,1 ,1 1 .40 
AFFX-LysX-5_at,6,H6at,2,60.70 
AFFX-LysX-5_at,6,Co!d,1 ,1 1 .90 
AFFX-LysX-5_at, 6, Cold, 2,9.70 
AFFX-LysX-M_at,1. Control, 1,75. 30 
AFFX-LysX-M_at.  1  ,Control,2,61 .90 
AFFX-LysX-M_at,  1  ,Heat,1 ,53.70 
AFFX-LysX-M_at.1  ,Heat,2,27.40 
AFFX-LysX-M_al,1 , Cold.  1 ,54.90 
AFFX-LysX-M_at,1  .Cold, 2,59.50 
AFFX-LysX-M_at,2,Control,1, 78.70 
AFFX-LysX-M_at,2,Control,2, 33.30 
AFFX-LysX-M_at,2,Heat,1,24.30 
AFFX-LysX-M_at,2,Heat.2.25.80 
AFFX-LysX-M_at,2,Cold,1 ,23.40 
AFFX-LysX-M_at,2,Cold,2,9.90 
AFFX-LysX-M_at,3,Control,  1 ,42.30 
AFFX-LysX-M_at,3,Control, 2,80.70 
AFFX-LysX-M_at,3,Heat.1 ,30.80 
AFFX-LysX-M_at,3, Heat, 2, 98.10 
AFFX-LysX-M_at,3,Cold,1, 20.30 
AFFX-LysX-M_at,3, Cold, 2,19.70 
AFFX-LysX-M_at,4,Control, 1,19.80 
AFFX-LysX-M_at,4,Control,2,24.20 
AFFX-LysX-M_at,4,Heat,1 ,27.80 
AFFX-LysX-M_at,4,Heat,2,50.50 
AFFX-LysX-M_at,4,Cold,1 ,20.20 
AFFX-LysX-M_at,4,Co!d,2,23.10 
AFFX-LysX-M_at,5,Control,1, 23.90 
AFFX-LysX-M_at,5,Control, 2,35.00 
AFFX-LysX-M_at,5, Heat,  1,15. 70 
AFFX-LysX-M_at,5,Heat,2,41.70 
AFFX-LysX-M_at,5,Cold,1 ,44.60 
AFFX-LysX-M_at,5.Cold,2,26.30 
AFFX-LysX-M_at,6,Control,1,20.10 
AFFX-LysX-M_at, 6, Control, 2,43.00 
AFFX-LysX-M_at,6,Heat,1 ,33.70 
AFFX-LysX-M_at,6,Heat,2,144.40 
AFFX-LysX-M_at,6,Cold,1, 78.20 
AFFX-LysX-M_at,6,Cold,2,85.30 
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AFFX-LysX-3_at,1  .Control, 1 ,8.80 
AFFX-LysX-3__at,1, Control, 2,46. 30 
AFFX-LysX-3_at,1  ,Heat,1 ,27.1 0 
AFFX-LysX-3_at,1 , Heat, 2,38.90 
AFFX-LysX-3_at,1  ,Cold,1 ,6.30 
AFFX-LysX-3_at,1  ,Cold,2,61 .50 
AFFX-LysX-3_at, 2, Control, 1,50.70 
AFFX-LysX-3_at,2, Control, 2,9.70 
AFFX-LysX-3_at,2,Heat,1 ,22.70 
AFFX-LysX-3_at, 2, Heat, 2,8.50 
AFFX-LysX-3_at,2,Cold,1,4.10 
AFFX-LysX-3_at,2,Cold,2,4.70 
AFFX-LysX-3_at,3,Control,  1 ,2. 1 0 
AFFX-LysX-3_at,3,Control,2,7.50 
AFFX-LysX-3_at,3,Heat,1 ,36.70 
AFFX-LysX-3_at,3,Heat,2,8.30 
AFFX-LysX-3_at,3,Cold,1 ,1 6.60 
AFFX-LysX-3_at,3,Cold,2,27.40 
AFFX-LysX-3_at.4, Control,!, 12.20 
AFFX-LysX-3_at,4,Control,2,10.70 
AFFX-LysX-3_at,4,Heat,  1 ,7.20 
AFFX-LysX-3_at,4,Heat,2, 56.80 
AFFX-LysX-3_at,4,Cold,1 ,44.40 
AFFX-LysX-3_at,4,Cold,2,30.50 
AFFX-LysX-3_at,5,Control,1, 19.90 
AFFX-LysX-3_at,5,Control,2,6.30 
AFFX-LysX-3_at,5,Heat,1,2.50 
AFFX-LysX-3_at,5,Heat,2,15.40 
AFFX-LysX-3_at,5,Cold,1 ,4.80 
AFFX-LysX-3_at,5,Cold, 2,27.10 
AFFX-LysX-3_at,6,Control,1, 36.00 
AFFX-LysX-3_at,6,Control,2,1 1 .60 
AFFX-LysX-3_at,6,Heat,1, 20.60 
AFFX-LysX-3_at,6,Heat,2,39.20 
AFFX-LysX-3_at,6,Cold,1 ,28.00 
AFFX-LysX-3_at,6,Cold, 2,8.00 
AFFX-PheX-5_at,1, Control, 1,8. 10 
AFFX-PheX-5_at,  1  .Control, 2, 5. 1 0 
AFFX-PheX-5_at,1  ,Heat,1 ,9.70 
AFFX-PheX-5_at,1, Heat, 2,6.50 
AFFX-PheX-5_at.1,Cold,1,5.10 
AFFX-PheX-5_at,1, Cold, 2,13.30 
AFFX-PheX-5_at,2,Control, 1,9.10 
AFFX-PheX-5_at,2,Control,2,10.50 
AFFX-PheX-5_at,2,Heat,1 ,1 2.50 
AFFX-PheX-5_at,2,Heat,2,10.70 
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AFFX-PheX-5_at, 2, Cold, 1,5.70 
AFFX-PheX-5_at,2,Cold,2.5.90 
AFFX-PheX-5_at,3,Control, 1,8.80 
AFFX-PheX-5_at,3,Control,2,1 1 .30 
AFFX-PheX-5_at,3,Heat,1 ,6.1 0 
AFFX-PheX-5_at,3, Heat, 2,6.40 
AFFX-PheX-5_at,3,Cold,1 ,1 1 .50 
AFFX-PheX-5_at,3,Cold,2,6.20 
AFFX-PheX-5_at,4,Control,1 ,4.90 
AFFX-PheX-5_at,4,Contro!, 2,3.40 
AFFX-PheX-5_at,4,Heat,1 ,6.90 
AFFX-PheX-5_at,4, Heat, 2,9.00 
AFFX-PheX-5_at,4.Cold, 1,7.20 
AFFX-PheX-5_at,4,Cold,2,7.80 
AFFX-PheX-5_at,5,Control,1 ,5.00 
AFFX-PheX-5_at,5,Control, 2,7.30 
AFFX-PheX-5_at,5,Heat,1 ,6.60 
AFFX-PheX-5_at,5,Heat, 2,5.10 
AFFX-PheX-5_at,5,Cold, 1,9.00 
AFFX-PheX-5_at,5,Cold,2,6.20 
AFFX-PheX-5_at, 6, Control, 1 ,7.40 
AFFX-PheX-5_at,6,Control,2,5.10 
AFFX-PheX-5_at,6,Heat,1 ,7.30 
AFFX-PheX-5_at.6,Heat,2,9.40 
AFFX-PheX-5_at,6,Cold,1 ,7.70 
AFFX-PheX-5_at,6,Cold.2,1 3.20 
AFFX-PheX-M_at,1  ,Control,1 ,8.20 
AFFX-PheX-M_at,1  ,Control,2,1 1 .60 
AFFX-PheX-W_at.1  .Heat,1 ,8.80 
AFFX-PheX-W_at,1.Heat,2, 12.30 
AFFX-PheX-M_at,  1 ,  Cold  ,1,12.60 
AFFX-PheX-M_at.1  ,Cold,2,1 6.80 
AFFX-PheX-M_at,2,Control,  1 , 1 7.90 
AFFX-PheX-M_at,2,Control,2,8.70 
AFFX-PheX-M_at,2,Heat, 1,5.90 
AFFX-PheX-M_at,2, Heat, 2, 16.20 
AFFX-PheX-M_at,2,Cold,  1 ,23.50 
AFFX-PheX-M_at,2,Cold.2,12.10 
AFFX-PheX-W_at,3,Control,1 ,8.90 
AFFX-PheX-M_at,3,Control,2,5.80 
AFFX-PheX-M_at, 3, Heat, 1,8.60 
AFFX-PheX-M_at,3,Heat,2,10.10 
AFFX-PheX-M_at,3,Cold,1 .4.30 
AFFX-PheX-M_at,3,Cold,2,8.10 
AFFX-PheX-M_at,4,Control,  1 , 1 0.70 
AFFX-PheX-M_at,4,Control,2,6.70 
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AFFX-PheX-M_at,4.Heat,1 ,9.80 
AFFX-PheX-M_at,4,Heat,2,9. 1 0 
AFFX-PheX-M_at,4,Cold,1 ,8.40 
AFFX-PheX-M_at,4,Cold,2,4.60 
AFFX-PheX-M_at, 5, Control,  1 ,4.90 
AFFX-PheX-M_at,5,Control,2,8.40 
AFFX-PheX-M_at,5,Heat,1 ,6.00 
AFFX-PheX-M_at,5,Heat, 2,5.90 
AFFX-PheX-M_at,5,Cold,1 ,5.20 
AFFX-PheX-M_at,5,Cold,2,6.90 
AFFX-PheX-M_at,6,Control,1 ,7.60 
AFFX-PheX-M_at,6,Control,2,7.70 
AFFX-PheX-M_at,6,Heat,1 ,1 7.90 
AFFX-PheX-M_at,6,Heat,2,28.20 
AFFX-PheX-M_at,6,Cold,1 ,27.00 
AFFX-PheX-M_at,6, Cold, 2, 13.10 
AFFX-PheX-3_at,1  .Control, 1 ,1 1 3.50 
AFFX-PheX-3_at,1, Control, 2, 28.80 
AFFX-PheX-3_at,1  ,Heat,1 ,65.70 
AFFX-PheX-3_at,1  ,Heat,2,1 50.70 
AFFX-PheX-3_at,1  .Cold,  1 ,73.20 
AFFX-PheX-3_at,1  ,Cold,2, 56.90 
AFFX-PheX-3_at,2,Control.  1,43.80 
AFFX-PheX-3_at,2,Control,2, 88.80 
AFFX-PheX-3_at,2,Heat,1 ,36.50 
AFFX-PheX-3_at,2,Heat,2,82.80 
AFFX-PheX-3_at,2,Cold,1, 98.90 
AFFX-PheX-3_at,2,Cold,2, 36.20 
AFFX-PheX-3_at,3,Control,  1 ,69.80 
AFFX-PheX-3_at,3,Control, 2,89.90 
AFFX-PheX-3_at,3,Heat,1, 95.60 
AFFX-PheX-3_at,3,Heat,2,61.90 
AFFX-PheX-3_at,3, Cold, 1,36.60 
AFFX-PheX-3_at,3,Cold,2,40.40 
AFFX-PheX-3_at,4,Contro!,  1 ,68.90 
AFFX-PheX-3_at,4,Controi,2,37.00 
AFFX-PheX-3_at,4,Heat,1 ,46.90 
AFFX-PheX-3_at,4,Heat,2, 109.90 
AFFX-PheX-3_at,4,Cold,1 ,58.40 
AFFX-PheX-3_at,4,Cold,2,91 .50 
AFFX-PheX-3_at,5,Control,1 ,1 07.80 
AFFX-PheX-3_at,5,Contro!,2,32.50 
AFFX-PheX-3_at,5,Heat,1,52.70 
AFFX-PheX-3_at,5,Heat,2, 37.60 
AFFX-PheX-3_at,5,Cold,1 ,1 02.90 
AFFX-PheX-3_at,5,Cold,2,1 1 8.40 
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AFFX-PheX-3_at, 6, Control,  1,85.70 
AFFX-PheX-3_at, 6, Control, 2, 99.50 
AFFX-PheX-3_at,6,Heat, 1,78.40 
AFFX-PheX-3_at,6,Heat,2,241.00 
AFFX-PheX-3_at,6,Cold,1 ,75.1 0 
AFFX-PheX-3_at,6,Cold,2,79.70 
AFFX-ThrX-5_at,  1  ,Control,1 ,82.30 
AFFX-ThrX-5_at.1, Control, 2,19.30 
AFFX-ThrX-5_at,1 , Heat, 1 ,23.70 
AFFX-ThrX-5_at,1  .Heat, 2, 26.40 
AFFX-ThrX-5_at,1  ,Cold,1 ,42.90 
AFFX-ThrX-5_at,1  ,Cold,2, 22.90 
AFFX-ThrX-5_at,2,Control,  1 ,26.20 
AFFX-ThrX-5_at,2,Control,2,34.30 
AFFX-ThrX-5_at,2,Heat,1 ,29.00 
AFFX-ThrX-5_at,2,Heat,2,31 .00 
AFFX-ThrX-5_at, 2, Cold, 1,30.40 
AFFX-ThrX-5_at,2,Cold,2,19.30 
AFFX-ThrX-5_at,3, Control, 1 ,1 7.30 
AFFX-ThrX-5_at,3,Control,2,12.00 
AFFX-ThrX-5_at,3,Heat,1 ,1 6.60 
AFFX-ThrX-5_at,3, Heat, 2, 18.80 
AFFX-ThrX-5_at,3,Cold,1 ,1 3.80 
AFFX-ThrX-5_at,3,Cold,2,17.70 
AFFX-ThrX-5_at, 4, Control,  1 ,35.60 
AFFX-ThrX-5_at,4, Control, 2, 17.00 
AFFX-ThrX-5_at,4,  Heat,  1 , 1 2 .50 
AFFX-ThrX-5_at,4,Heat,2,18.70 
AFFX-ThrX-5_at.4,Cold.1, 20.40 
AFFX-ThrX-5_at,4,Cold.2,9.60 
AFFX-ThrX-5_at,5,Control,1 ,21 .80 
AFFX-ThrX-5_at,5,Control,2,15.00 
AFFX-ThrX-5_at,5,Heat,1,8.00 
AFFX-ThrX-5_at,5,Heat,2. 18.30 
AFFX-ThrX-5_at,5,Cold, 1.1 3.90 
AFFX-ThrX-5_at,5,Cold, 2,15.40 
AFFX-ThrX-5_at,6,Control, 1,14.60 
AFFX-ThrX-5_at,6, Control, 2,13.60 
AFFX-ThrX-5_at, 6, Heat, 1,17.50 
AFFX-ThrX-5_at,6,Heat,2,17.50 
AFFX-ThrX-5_at,6,Cold,1 ,21 .30 
AFFX-ThrX-5_at,6,Cold,2,17.70 
AFFX-ThrX-M_at,1 , Control,  1 ,5.90 
AFFX-ThrX-M_at,1,Control,2, 56.70 
AFFX-ThrX-M_at,1  ,Heat,1 ,62.30 
AFFX-ThrX-M_at,1  ,Heat,2,59.00 
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AFFX-ThrX-M_at,1,Cold.1, 54.10 
AFFX-ThrX-M_at,  1 , Cold, 2, 53.20 
AFFX-ThrX-M_at, 2, Control, 1,76.90 
AFFX-ThrX-M_at,2, Control, 2,38.70 
AFFX-ThrX-M_at,2,Heat,  1 , 1 7.20 
AFFX-ThrX-M_at,2,Heat,2,78.60 
AFFX-ThrX-M_at,2,Cold,1 ,62.90 
AFFX-ThrX-M_at, 2, Cold, 2,43.40 
AFFX-ThrX-M_at, 3, Control, 1,52.00 
AFFX-ThrX-M_at,3,Control,2,60.20 
AFFX-ThrX-M_at,3,Heat,1, 74.90 
AFFX-ThrX-M_at,3,Heat,2,47.20 
AFFX-ThrX-M_at,3,Cold.  1 , 1 3.00 
AFFX-ThrX-M_at,3,Cold, 2,20.60 
AFFX-ThrX-M_at,4,Control, 1,66.70 
AFFX-ThrX-M_at,4,Control,2,6.10 
AFFX-ThrX-M_at,4,Heat,1 ,1 1 .20 
AFFX-ThrX-M_at,4,Heat,2,18.10 
AFFX-ThrX-M_at,4,Cold,1 ,53.20 
AFFX-ThrX-M_at.4.Cold,2,29.70 
AFFX-ThrX-M_at,5,Control,1, 14.50 
AFFX-ThrX-M_at,5,Control, 2,42.80 
AFFX-ThrX-M_at,5,Heat,  1 ,44.90 
AFFX-ThrX-M_at,5,Heat,2, 60.70 
AFFX-ThrX-M_at,5,Cold,1 ,40.30 
AFFX-ThrX-M_at,5,Cold,2,34.80 
AFFX-ThrX-M_at,6,Control,1, 60.80 
AFFX-ThrX-M_at,6,Control,2,17.80 
AFFX-ThrX-M_at,6,Heat,1 ,68.80 
AFFX-ThrX-M_at,6,Heat,2,38.20 
AFFX-ThrX-M_at,6,Cold,  1 ,55.30 
AFFX-ThrX-M_at,6,Cold,2,53.10 
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TEST  TWO-WAY  ANOVA  OUTPUT  FILE 


Probe 

p  value  for 
condition 

p  value  for  time 

p  value  for 
interaction 

AFFX-BioB-5_at 

0.05794984 

0.13501568 

0.24353257 

AFFX-BioB-M_at 

0.35676659 

0.30461054 

0.59663903 

AFFX-BioB-3_at 

0.06939791 

0.21827141 

0,41437524 

AFFX-BioC-5_at 

0.24299456 

0.30101028 

0.52575993 

AFFX-BioC-3_at 

0.22649861 

0.31565169 

0.41011583 

AFFX-BioDn-5_at 

0.11384801 

0.19242631 

0.28908168 

AFFX-BioDn-3_at 

0.16022979 

0.16946863 

0.27113168 

AFFX-CreX-5_at 

0.13472166 

0.09849027 

0.43169993 

AFFX-CreX-3_at 

0.13544086 

0.08311014 

0.42987069 

AFFX-DapX-5_at 

0.95629408 

0.82901326 

0.71919154 

AFFX-DapX-M_at 

0.94023323 

0.68961177 

0.08185544 

AFFX-DapX-3_at 

0.30338196 

0.67242893 

0.10017086 

AFFX-LysX-5_at 

0.84402251 

0.77746585 

0.86836897 

AFFX-LysX-M_at 

0.79542281 

0.17446481 

0.16022286 

AFFX-LysX-3_at 

0.70302937 

0.56640906 

0.55522731 

AFFX-PheX-5_at 

0.63177669 

0.80000401 

0.75977661 

AFFX-PheX-M_at 

0.57126931 

0.8818067 

0.12279424 

AFFX-PheX-3_at 

0.53834261 

0.46861711 

0.0983472 

AFFX-ThrX-5_at 

0.52682431 

0.14079342 

0.16039673 

AFFX-ThrX-M_at 

0.70419305 

0.59522396 

0.76935162 

No.  of  Genes  (P  <=  0.01 ) 

0 

0 

0 
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